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FOREWORD 


This  report  was  prepared  by  the  Lockheed-California  Company  under  contract 
DOT-FA75-WA-3707 .  The  report  contains  a  partial  description  of  the  effort 
performed  as  Task  II  of  a  three  task  effort.  Task  II  covers  the  period  from 
January  1978  to  November  1978.  The  work  was  administered  under  the  direction 
of  the  Federal  Aviation  Administration  with  H.  Spicer  as  Technical  Monitor. 

The  project  leader  was  Oil  Wittlin,  of  the  Lockheed-California  Company. 
Important  contributions  were  made  to  the  program  by  the  Cessna  Aircraft 
Company  under  the  direction  of  D.J.  Ahrens  and  W.B.  Bloedel.  W.L.  LaBarge, 
of  the  Lockheed-California  Company,  prepared  the  Programmer's  Manual  with 
the  assistance  of  0.  Wittlin,  M.S.  Camon,  H.  Weinberger  and  D.  Seiki. 

The  Lockheed  effort  was  performed  under  the  supervision  of  J.E.  Wignot 
(Dynamic  Loads  Group). 
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SUMMARY 


This  document  contains  a  description  of  program  KRASH  as  modified  under  con¬ 
tract  DOT- FA75-WA- 3707.  Included  in  the  Programmer's  Manual  are  the  following 

sections : 

•  Program  KRASH  System  Requirements 

•  Input  Data  Deck 

•  Demonstration  Problem 

This  Programmer’s  Manual  is  one  of  a  series  of  operational  documents  for 
program  KRASH.  In  this  manual,  pertinent  information  is  supplied  which  will 
facilitate  bringing  the  program  to  an  operational  status  on  the  user's 
computer  system.  This  manual  has  been  established  in  such  a  manner  that  it 
can  be  updated  as  more  data  becomes  available.  The  subject  material  contained 
within  each  section  can  be  expanded  or  revised,  as  necessary,  without  affecting 
the  other  sections.  Each  section  contains  its  own  numbering  system  which 
facilitates  the  task  of  updating  the  document. 
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SECTION  1 


INTRODUCTION 


Program  KRASH  is  a  digital  computer  program  which  predicts  the  struc¬ 
tural  response  of  vehicles  to  multidirectional  crash  environments.  The 
program  computes  the  time  histories  of  N  interconnected  masses,  each  allowed 
six  degrees  of  freedom,  defined  by  internal  coordinates  X^,  Y.,  and 
Eulerian  angles  ,  e^,  X  ,  when  i  =  1,2,...N.  Euler's  equations  of  motion 
are  written  for  each  mass.  The  equations  of  motion  are  integrated  numerically 
to  obtain  velocities,  displacements,  and  rotations.  Gravitational  forces, 
internal  forces  and  moments,  and  external  forces  are  computed.  For  small 
deflections,  a  linear  analysis  is  obtained,  and  for  large  deflections,  general 
plastic  deformation  is  allowed.  The  program  provides  for  unloading  and  subse¬ 
quent  reloading  along  a  linear  elastic  line. 

Program  KRASH  describes  the  interaction  between  a  series  of  massless 
interconnecting  structural  elements  and  concentrated  rigid  body  masses  to 
which  the  structural  elements  are  attached  at  their  ends  with  the  appropriate 
end  fixity  (pinned,  fixed).  The  structural  elements  can  be  connected  between 
node  points  which  are  offset  from,  and  rigidly  attached  to,  selected  mass 
points.  The  interconnecting  elements  represent  the  stiffness  characteristics 
of  the  structure  between  the  masses.  The  masses  can  translate  and  rotate  in 
all  directions  under  the  influence  of  the  external  forces  (i.e.,  gravity, 
aerodynamic,  impact),  as  well  as  the  constraining  internal  element  forces. 

The  movement  of  the  masses  results  in  changes  in  the  relative  distortion  of 
the  structural  elements  and,  in  turn,  results  in  a  new  set  of  element  forces 
acting  throughout  the  system. 
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Computer  Program  KRASH  has  the  capability  to: 

•  Define  the  response  of  six  degrees  of  freedom  (DOF)  at  each  repre¬ 
sentative  location,  including  three  translations  and  three  rotations. 

•  Determine  mass  accelerations,  velocities,  and  displacements  and 
internal  member  loads  and  deformations  at  each  time  interval. 

•  Provide  for  general  nonlinear  stiffness  properties  in  the  plastic 
regime,  including  different  types  of  load-limiting  devices,  and 
determine  the  amount  of  permanent  deformation. 

•  Define  how  and  when  rupture  of  an  element  takes  place  and  redistribut 
the  loading  over  the  structural  elements  involved. 

•  Define  mass  penetration  into  an  occupiable  volume. 

•  Define  the  volume  change  due  to  structural  deformations  of  an  occu¬ 
piable  volume. 

•  Provide  for  ground  contact  by  external  structure  including  sliding 
friction  and  a  nonrigid  ground  surface. 

•  Include  internal  structural  damping. 

•  Include  a  measure  of  injury  potential  to  the  occupants;  for  instance, 
the  probability  of  spinal  injury  indicated  by  the  Dynamic  Response 
Index  (DRI). 

•  Determine  the  distribution  of  kinetic  and  potential  energy  by  mass 
item,  the  distribution  of  strain  and  damping  energy  by  beam  element, 
and  the  crushing  and  sliding  friction  energy  associated  with  each 
external  spring. 

•  Determine  the  vehicle  response  to  an  initial  condition  that  includes 
linear  and  angular  velocity  about  three  axes  and  any  arbitrary  vehicl 
attitude  and  position. 

•  Provide  a  measure  of  the  airplane  eg  velocity  by  means  of  transla¬ 
tional  momentum  relationships. 

•  Analyze  an  impact  into  a  horizontal  ground  and/o-  an  inclined  slope. 

•  Provide  a  measure  of  the  internal  stress  state  of  internal  beam 
elements. 

•  Analyze  a  mathematical  model  containing  up  to  80  masses  and  150 
Internal  beam  elements. 

•  Treat  up  to  180  nonlinear  element  degrees  of  freedom. 
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A  comprehensive  discussion  of  the  theoretical  development  of  program 
KRASH  is  presented  in  Reference  (1).  A  discussion  of  program  KRASH  input- 
output  formats  as  well  as  modeling  techniques  and  applications  are  included 
in  Reference  (2).  Reference  (3)  provides  a  discussion  of  design  guidelines 
which  can  be  used  in  conjunction  with  program  KRASH  in  the  structural  crash- 
worthiness  analysis  of  general  aviation  airplanes. 

The  last  in  the  series  of  operational  documentation  for  program  KRASH  is 
this  Programmer's  Manual.  In  this  manual,  pertinent  information  is  supplied 
which  will  facilitate  bringing  the  program  to  an  operational  status  on  the 
user's  computer  system.  The  manual  contains  a  discussion  of  the  program's 
system  requirements,  functional  organization,  and  a  demonstration  problem 
which  can  be  used  to  verify  the  operation  of  the  program. 
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SECTION  2 


PROGRAM  KRASH  SYSTEM  REQUIREMENTS 


The  basic  control  card  deck  consists  of  the  necessary  job  control  cards, 
the  program  KRASH  source  cards,  and  the  data  input  card  deck.  Here,  "card 
deck"  and  "cards"  are  to  be  interpreted  in  their  broadest  sense.  They  include 
any  medium  by  which  instructions  are  transmitted  to  the  computing  system,  e.g., 
physical  cards,  tape,  remote  terminal  instructions,  etc. 

2.1  JOB  CONTROL  DECK 

Due  to  the  wide  variety  of  computer  systems  available  and  individual  com¬ 
puter  installation  requirements,  a  complete  and  detailed  listing  of  all 
possible  job  control  language  (JCL)  cards  is  not  feasible.  Format  details  of 
the  JCL  instructions  will  have  to  satisfy  the  requirements  of  the  specific 
computing  system  in  use.  Irrespective  of  the  details,  the  JCL  cards  must 
provide  for  the  following  program  KRASH  requirements: 

•  Using  the  overlay  structure  illustrated  in  Figure  2-1,  a  core  space 

equivalent  to  600K  bytes  of  IBM  360  core.  (Note  that  the  IBM  360 

is  a  32-bit  word  machine  with  4  bytes  per  word). 

•  The  temporary  logical  data  storage  units  discussed  in  Section  2.2. 

•  A  permanent  data  storage  unit.  See  Section  2.2. 

•  The  compile  and  link  instructions  necessary  to  execute  a  FORTRAN  IV 
source  deck.  For  the  examples  shown  in  this  manual,  a  FORTRAN  H 
compiler  (level  21.8)  with  an  optimization  level  of  2  (OPT  22)  has 
been  used . 

2.2  LOGICAL  UNITS 

Calls  to  the  eleven  logical  storage  units  listed  in  Table  2-1  can  be 

made  from  program  KRASH.  Of  the  units  listed,  all  but  units  5,  6,  and  21  are 
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Figure  2-1.  Overlay  Map 


TABLE  2-1.  LOGICAL  UNIT  SUMMARY 


considered  temporary  data  storage  devices.  Of  the  temporary  units,  unit  10 
is  required  for  each  run.  The  remaining  temporary  units  are  required  only 
when  summary  plots  are  requested.  See  Section  2.1  of  Reference  (2).  Units 
5  and  6  are  standard  stored  input  and  output  units.  Unit  21  is  a  permanent 
data  storage  unit  and  is  called  whenever  the  restart  option  is  exercised 
(see  Section  3.3). 

2.3  INPUT  DATA  DECK 

The  input  data  required  to  run  program  KRASH  consists  of  a  series  of 
cards  which  control  the  operation  of  the  program,  set  the  initial  conditions 
of  the  analyses,  and  define  the  structural  characteristics  of  the  model.  A 
complete  description  of  the  cards  which  make  up  this  deck  is  presented  in 
Section  2  of  Reference  (2).  All  of  the  data  supplied  on  the  cards  of  this 
deck  are  read  from  unit  5. 

Special  attention  is  directed  to  cards  3  and  6  of  the  input  data  deck 
of  program  KRASH.  See  Figure  2-1  of  reference  (2).  The  data  supplied  on  these 
cards  are  used  when  the  restart  option  is  exercised.  The  restart  option 
requires  a  check  point  run  to  be  made  during  which  state  vector  data  are  saved 
on  unit  21  at  specified  times.  In  subsequent  runs  the  data  can  be  retrieved 
from  unit  21  and  the  program  restarted  at  any  one  of  the  times  used  in  the 
check  point  run.  A  description  of  the  data  management  details  related  to  the 
storage  and  retrieval  of  these  data  is  presented  in  Section  3.3. 

2.4  RUNNING  PROGRAM  KRASH 


2.4.1  Program  KRASH  Execution 

Figure  2-2  presents  a  copy  of  the  program  deck  used  at  the  Lockheed 
California  Company  to  compile,  link,  and  execute  program  KRASH.  The  program 
deck  assumes  that  program  KRASH  source  is  on  an  unlabeled  9-track  tape  with 
card  image  records  (i.e.  RECFM  -  F,  LRECL  ■  80,  BLKSIZE  -  80);  note  that  this 
is  the  standard  format  of  the  KRASH  source  tape  as  released  by  the  Lockheed 
California  Con^any.  If  the  program  KRASH  source  is  on  cards,  then  changes 
to  the  exanple  program  deck  will  be  necessary  to  reflect  this  fact. 
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2.U.2  Restart  Tape  Initiation 


Before  the  restart  option  of  program  KRASH  tan  be  exercised,  a  save 
tap,  must  be  reserved  and  Initialised  as  described  In  Section  3.3.  The  pro- 
gre.  deck  ah™  In  Figure  2-3  he,  been  os.d  at  the  Lodthe.d  California  Compare 

to  Initialise  a  re.t.rt  tape.  The  tap.  In  the  exaeple  Is  an  IBM  Standard 
Label  tape. 


//KRA.VH  JOB  NNNN  ,  NAMI  ,  RE  G  I  ON  *600K  ,  1  I  HI  I  S 
/•  /('  EXE C  ROM  IEKAAOO  ,  F  ARM  ■  •  MAP  B<  D  NUDE  CK  ,  OP  I  2  1 D  ’ 

//SY SPRINT  DP  SYSOIJT -A 
/SYSPUNCH  DI)  SYSOUT=B 

/ / S Y S LIN  D I)  I) S  NAME  *«.COMP, SPA C E  =  < 3 2 0 0  ,  (  1  0 O ,  5 0  )  >  ,  U N  I  I  3 3 3 0  , 
//  DCB=  (  RECFM  FB  ,  LRECL  =£JO  ,  BL.KSI/E  -  3200 )  ,  DISC  (NEW  .  PAS’S*  > 
//SYSUT2  DD  UN  I  I  -  3330  ,  SPAC  E=  <  TRK  ,  (5,3)  ) 

//SYSABEND  DD  SYSOUT  A  ,  DCB=  (  BLKSI  ZF  =800  > 

//**  ¥  ********  *  *  *  *  «*  *  *  *.**¥  *  *.  **  #*  *  *  *  *  r  *  *  ***  ,  *;  >•*  »  *  >  * 

//# 

//*  THE  KRASH  SOURCE  TAPE  AS  RELEASED  IS  AN  UNI  AB!  I  I  T)  *. 
■'■■*  V  TRACK  TAPE  WITH  CARD  IMAGE  RECORD.". 

•’  /  * 

/  /  *  K  *  *******  *  *  *  *  K  *  *  ■  *  *  *  *  *  *  *  *  *  *  »  *>'  *  ft  *  >  ft  *  K  *  >  >  K  «•  *t>  *  *»  *  »»<•»>'»  »  » 

//SYS  IN  DD  UNI  F=240M,I.  ABEL  (  ,  NL  )  ,  Df  B  v  RE  Cl  M  i  .1  II  Ci  HO,  H*  I 
V  D ISP -(OLD, KEEP) 

//I.  EXEC  PGM  TEWI  .  I  'ARM  ••  1  L  1ST  ,  MAP ,  OVI  Y  ’  .  COND ( '■  I  I  ,  C  >  ,  T  f  Ml 
/  /  S  Y  ,V  I  ■  R I N  T  D  D  S  Y  S  0  (J  T  -  A 

/ / r 0 Ml  I  LF  DD  D S N A M £ ~ A 0 0 M P  D  I S P  =  (  01.  T> ,  D L L •  FT  ■ 

/■'SYSUT 1  DD  DCB-  (  BL  KSIZE=7294  >  ,  ST  A  Cl  (7274,  (40. 20  >  >  ,  UNI  T 
/  /  S  Y  S  L  I  B  D  D  D  S  NAME  S  Y  S  I  .  F  0  R  1 1  I  B  ,  D  J  S  P = S  H  P 
//  SYSI.  MOD  DO  DSN  AMI:  =  &  FEME'  ( GOMODUL  E  )  ,  DISP=<NEW,  I  ASS  >  , 

//  UNI  T  =3330  ,  SPACE  (  7294  .  ( 40 , 20 . 1  )  )  ,  DCB  (  JHI.KSI  7E=  72V 4  > 
■••SYSLIN  DD  k 
I  NT 'I  ODE  COMP  1 1  F 
ENTRY  MAIN 
OVERl AY  A 

I  NSERT  OF  R  I.V,  ACCEL T  ,  INTTRP.CF  01 'Cl 
OVERLAY  A 

l  ysi:  r;  I  HI  nmod  ,  INPUT  ,  FNI'RNT  ,  ECHO 
OVERLAY  A 

!  N  S  i:  R  T  R  S  0  U  T  ,  D  A  T  I  N  ,  R  S I N  .  D  A  T  0  U  T 
I  NS I.  R  F  JOBAC  ,  ('DATE 
I  NSF  R  I  RC  ,  DOA  I .  I ,  I C  .  MA  T  Ml  l|._  ,  PR  I  N  T 
I NSt  R  T  PAPL.T  ,  V SHELL  ,  SHI.  LI.  M  ,  SHELLY. 

OVERLAY  A 

INSERT  F  RF  PI  I  ,  PI  OFT 

//G  EXEC  PGM  *>.  .1,  „  SYSLMOD  ,  COND=  (  (5,LT,C>  ,  (5  .  IT  .1  )  )  ,  REGION 
//F  T  01  r- 001  DD  DCB=  (RECFMVB  BEK SIZE =000, 1  PI  Cl.  =776), 

//  UNI  F =3330  .  SPACE  =  <  I  RK  ,  ( 1  0  ,  5  >  > 

//F  T02E001  DD  DCB  -  (RE  Cl  M  VB  ,  BL  R  SIZES"  000  ,  l.  RE  CL  =  796)  , 

//  UNI T=3330 , SPACE* ( TRK , < I  0 ,5)) 

//I  T  031  001  DD  DCB  (RECFM  VB  ,  BLKSIZE  =000  ,  LRI  CL  v°7,)  , 

//  UN I T= 3330 , SPACE* ( TRK , (10,5) ) 

//FT04I  001  DD  DCB  :  (  RECF  M-  VB  ,  BLKSI ZE  =  800  ,  ERE  Cl.  =  796  )  , 

//  UNI T =3330 , SPACE = ( TRK , (10,5)) 

//FI OAF 001  DD  SYSOUT=A 
//FT07F001  DD  SYS0UT=B 

//FT08F001  DD  DCB=  ( RECF  M=VB  ,  BEK  SIZE  =  800  ,  L  RECL  =  79Z, )  , 

//  UNIT=3330 , SPACE* <  TRK ,(10,5)) 

Figure  2-2.  Example  Compile,  Link,  and  Execute 
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a  •  1  V  4  i  C  V1 


.',TTC?F00i  DP  PCD  ■'  Rl  (7  H  VP,  P!  i  I  71  POO  ,  I.  KH.I.  7V7  >  , 

//  UN  1 1  =  3330  ,  SPACt  -  <  T  RK  ,  ( 1  0  ,  M  > 

//FT  1  or- 001  DD  DCB=<RECFM=VB,BU  SIZE  TTOO ,  l  RE  CL  776  )  . 

/./  UN  I T  =  3330 , SPACE* <  TRK  ,  <  1  O  ,*i\\ 

//FT  1  1  FOOT  DD  DCB*  <Rfc.CFM=VB ,  Bl  I  S  TZf  B0>>  ,  l  R!  IT  ~~7'>’6  >  , 

//  UNT  1=3330  ,  SPACE  <  TRK  ,  00,5  >  > 

//FT  1  2F001  DD  DCB~  < RECFM=VB ,  BLK SIZE *800 ,  L.RECL ~7 '*6  >  . 

//  UNIT-3330 , SPACE* < T RK , <10, 5) ) 

//FT1  3F001  DD  BCB*  (  RECFM=VB ,  Bl  I  SIZE*SOO  .  I  RFC!  =79/. ) 

//  UN  .1  T= 3330  ,  SPACE*  <  TRK  ,  (  1  0 ,5)  ) 

/  /*#•***#■*****•  *  *  #  *;  >:  *  ■  *:  *  *■  *  *■*»•*■*  »*!•*>  **»>)>#*  *  X.  »  ■»  »  *  n  *  *  * ■  *  y  *  )>  «>  »  *  »  «  »  ii  (•  » 


//* 

//*  THf:  '\R  AS’M  CHECKPOINT/RE^T  ART  T  Ai  i  I  ,*  *  TO  Bl  H00» !  T I  o  (H 

//*  UNI  T  21.  DAT  A  SE  T  NANI  I  .v  1  M-  A/*M  .  DAT  A  '  .  TAPE  NMN.W  f.  f  ' 

//*  402512* 

/  /  H 


p 

* 

p 

M 


i  H  ¥  ¥:  ¥■  ¥  ¥  ¥  ¥;  ¥¥¥¥¥¥¥¥  *  )•  P  ¥ 


>  •"*»  >•  P  ¥  ¥:  P  P  H  H  -H  *  P  P  -P  ♦*  P  •»  ¥  P 


//FT21FO01  DD  UN  T  '1=2400,  DC»=<L  RECL  1  0  >4  .  T<!  '  :-l  71  ’  0:.’*  .  Rl  (  I  f)  VP' 

//  VOL*SER*40251 2 , LABEL* < ,SL> . D1 SP  0!  P  DSN=KRASH.DAIA 
//SYSUDUMP  DD  SYSOUT*A 
//TT05F001  DD  * 


< DATA  CARD  DECK ) 
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Drr- 


/  / 1 N I  I  TAPE  JOB  NNNN  ,  NAME  ,  REG  I  ON  ~  1  *01.  ,  T  I  Ml  2 


!•  *  *  ft  *  ******  *.  *.  *  V  >  * 


//*  THI S  l S  A  SAMPLE  JOB  OF  HOW  THI  KRASH  CHECKPOJ  N  T  /RE."  TART  * 

//*  TAPE  CAN  BE  INITIALIZED.  IBM  STANDARD  LABEL  ASSUMED  *. 

//*  NOT F  THAT  THE  TAPE  IS  A  V  TRACK  TAPE.  * 

/  /  *  * 

//<  EXEC  PGM* IEK AAOO  ,  REGION- 1  AOK  ,  FARM '•= '  MAP  ,  BCD  .  NO  DECK  ,  OPT- 1  ,  ID’  ,  TIMi.  I 
/•SYSPRINT  DD  SYSOUT-  A 
//SYSPUNCH  DD  SYSOUT ~B 

//SYS l  1 N  DD  DSNAME  ■  &COMP , SPACE  < 3200 . < 1 00 . MO ' 1 , UNIT =3330 , 

DCB* ( Rf  Cf  M  r B , LRECL =80 , BLKS I /E-  3200 > , D I SP= < NEW , PASS ) 

//SYSUT2  DD  UNIT-3330 (  SPACE* < TRK  ,  <3,  3)  1 
//SYSABE ND  DD  SYSOUT -A , DC B  < BLKSIZE=800 ) 

/SYS IN  DD  * 

ENI)  FILE  4 

S  I  HI  - 

END 

//I  E XEC  PGM  -  IEWI  ,  I  ARM  '  LIST  ,  MAP  ’  ,  FOND  <  5 , 1.  I  ,C>  ,  REGION  I  AOK,  I  I.MI  1 
//SYSPRINT  DD  SYSOUT  A 

//COMP.T  I.  I  DD  DSNAME  ’  SCOMI '  ,  DISC  <01  D  ,  D!  II  II) 


//SYSIM  I  DD  DCM*<  BLKSI  i'l  =  72V4  )  ,  SPACE  *  <  72  V  4  ,  <  40, 20  )  >  ,  UN  I.  P  3  -'30 

//SY.'LIB  DD  DSNAME  -SY  SI  EOI"  T I  IB  ,  DTSP~  SUP 

//:••»•  SI  MUD  DD  DSNAMI  &TL«P<  GOMODULE  )  ,  D.l  SP~  <  NEW,  PASS  >  . 

/ /  I  IN  I  I  3330  ,  SPACE*  <  729 4  ,  <  40 ,20,1  11,  DC  B  <  B!  PS T  /I  .•?'.* 4  ' 

••  SYSI  IN  DD  * 

INI  I  UDf  COMPILE 
f  MIR Y  MAIN 

G  I  XI  C  PGM  f  ,  I.  .SYSI.  MUD,  POND  <<M,  I..  I  ,  C)  ,  ,  LI  ,1...  1  >  ,  REGION  I  A  OP  ,  T  I  M1  1 

*  ******  *  *  ******  *  *  **  *  *  *:■  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *■**  *  V  ***********  n  ******  * 

/  /  *  * 

//*  MIL  KRASH  CHECKPOINT /RESTART  TAPE  IS  HI  BE  MOUNT  I  0  UN  » 

»  UNIT  4  DATA  SET  NAMI  IS  '  KRASH .  DA  I A '  ■  TAPI  NUMBER  I"  > 

>■>*  402312.  * 

•'  f  *  ¥ 

****************************  *  *  *  *  *  *  *  *  *************************** 

/I  I04E0O1  DD  DSN*  PR  ASH  .  DAT  A  ,  UNI  T= 2  400  ,  D1SP--  <  NE  W  ,  PI  I  pi, 

••/  DCB  <  RE  Cl  M  VB,  I.RECL  1  004  ,  BLKSIZE-3024  )  .  VOL  SEE  402’  >1  2  ,  L  ADI-  I  -<  .Si  > 
/ / 1  T  O A 1001  D  D  S Y SOU  T  A 
//SYSUDUMI  DD  SYSOUT=A 


Figure  2-3.  Example  Restart  Tape  Initialization 
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SECTION  3 


PROGRAM  KRASH  FUNCTIONAL  ORGANIZATION 


Program  KRASH  is  divided  into  three  basic  functional  blocks:  computa¬ 
tional  analysis  and  data  processing;  internal  data  storage  and  transfer; 
and  external  data  storage  and  retrieval.  Each  of  these  blocks  is  discussed 
in  the  following  sections  of  this  manual. 

3.1  COMPUTATIONAL  ANALYSIS  AND  DATA  PROCESSING 

The  computational  analysis  and  processing  of  data  is  accomplished  by 
program  KRASH  as  illustrated  in  the  flow  chart  shown  in  Figure  3-1.  The 
rectangular  blocks  represent  calls  to  functional  subroutines.  All  other 
blocks  represent  functional  operations  and  logical  branching  controlled 
within  the  MAIN  segment  of  the  program. 

3.1.1  Subroutine  Call  Cross  Reference 

Calls  to  various  subroutines  are  made  throughout  program  KRASH.  A  cross- 
reference  map  identifying  the  calling  routine  or  subroutine  and  the  called 
subroutine  is  shown  in  Table  3-1.  Included  in  the  map  are  the  IBM  360  core 
requirements  for  each  of  the  routines  listed.  Attention  is  drawn  to  the 
fact  that  these  size  requirements  are  applicable  only  when  the  H  version 
(OPT  ■  2)  of  the  FORTRAN  compiler  is  used.  The  use  of  other  compiler  ver¬ 
sions  or  options  will  result  in  different  size  requirements. 


— -—?-n 


TABLE  3-1.  SUBROUTINE  CROSS  REFERENCE  SUMMARY 


3-5 


RSIN 

RSOUT 

SHF.LLM 

SHELLX 


3.1.2  Subroutine  Functional  Description 


Described  in  this  section  are  the  functional  operations  of  each  subrou¬ 
tine  along  with  definitions  of  the  symbols  which  compose  the  subroutine  call 
statement  argument  list.  When  possible  these  symbols  are  referenced  to  the 
corresponding  symbols  used  in  the  theoretical  development  presented  in 
Reference  (1). 

3. 1.2.1  MAIN 

The  MAIN  segment  of  program  KRASH  serves  primarily  as  a  control  routine 
for  the  program.  Processing  of  data  is  controlled  by  calls  to  the  various 
functional  subroutines  included  in  the  program.  In  addition  to  the  above, 
the  MAIN  segment  contains  the  predictor-corrector  numerical  integration 
algorithm  used  by  program  KRASH.  The  predictor-corrector  method  used  is 
discussed  in  detail  in  Section  1.3.13  of  Reference  (1).  At  the  end  of  each 
run,  MAIN  controls  the  printing  of  various  summary  data  which  recap  the 
significant  events  of  the  run.  Examples  of  these  summaries  are  shown  in 
Figure  3.2. 

Comment  cards  are  included  in  the  MAIN  source  deck  which  refer  to 
RUNMOD  options  3  and  4.  These  options  are  not  discussed  in  the  User's  Manual, 
Reference  (2),  and  are  not  intended  for  production  use  at  this  time.  Option  3 
is  set  aside  for  future  use  as  a  model  data  card  punch  option.  Option  4  is 
an  option  which  has  been  set  up  to  aid  in  the  debugging  of  the  program  dur¬ 
ing  its  development. 

A  number  of  energy  error  checks  are  controlled  from  MAIN,  see  Section  3.4. 
When  any  of  these  checks  fail,  control  is  passed  to  the  summary  print  sec¬ 
tion  of  MAIN  for  execution  and  subsequent  run  termination. 
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Figure  3-2  Example  End  of  Run  Summary  Tables 


3. 1.2. 2  EULER 


Subroutine  EULER  generates  the  transformation  matrix  defined  by  equa¬ 
tion  1-2  of  Reference  (1).  The  angles  PHI,  THETA  and  PSI  are  used  to  form 
the  nine  elements  of  transformation  matrix  A.  The  elements  of  the  matrix 
are  returned  to  the  calling  routine  via  the  argument  list. 


CALLING  STATEMENT  ARGUMENT  LIST 

1 

Symbol 

Mnemonic 

Ref.  (1) 

Type 

A 

A 

R*8 

Transformation  matrix 

PHI 

4> 

R*8 

) 

THETA 

e 

R*8 

(  Euler  angles  defined  by  equa- 

PSI 

* 

R*8 

\  tion  1-2,  Reference  (1) 

3. 1.2. 3  MATVEC 

Subroutine  MATVEC  performs  the  vector  matrix  multiplication  A*V  if  ISW 
is  equal  to  zero,  otherwise  the  vector  matrix  multiplication  AT*V  is  per¬ 
formed.  The  product  of  the  multiplication  is  returned  to  the  calling  sub¬ 
routine  by  P. 


CALLING  STATEMENT  ARGUMENT  LIST 

Mnemonic 

Symbol 
Ref.  (1) 

Type 

A 

R*8 

Insert  matrix 

V 

R*8 

Input  matrix 

P 

R*8 

Resulting  output  matrix 

ISW 

1*4 

Transposition  flag 

A  blank  in  the  symbol  (Reference(l))  column  indicates  that  there 
is  no  equivalent  symbol  defined  in  the  KRASH  theory  which  can  be 
directly  related  to  the  particular  Mnemonic  term. 
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3. 1.2.4  ACCELT 

Subroutine  ACCELT  enters  the  acceleration  input  table  stored  in  the  common 
region  named  INAC.  Appropriate  levels  are  extracted  from  the  table  based  on 
the  time  parameter  TIME.  These  levels  are  passed  to  the  subroutine  INTERP  for 
interpolation.  The  interpolated  values  are  returned  to  the  calling  routine  by 
XI,  YC,  ZC,  A1 ,  A2,  and  A3.  The  counters  NPTSP  and  NPTS  are  set  equal  to  zero 
each  time  through  DERIV. 


CALLING  STATEMENT  ARGUMENT  LIST 

Mnemonic 

Symbol 
Ref.  (1) 

Type 

MASDEX 

1*4 

Mass  number 

l 

i 

XI 

u 

R*8 

> 

Time  derivative 

rxi 

YC 

V 

R*8 

>-  of  translational  - 

►  axis 

ZC 

• 

w 

R*8 

velocity  along 

A1 

p 

R*8 

I 

Time  derivative 

r x] 

A2 

q 

R*8 

►  of  rotational  -i 

Y 

►  axis 

A3 

• 

r 

R*8 

a 

velocity  about 

L  z  J 

TIME 

t 

R*8 

Integration  time  step 

NPTSP 

1*4 

counter,  previous 

NPTS 

1*4 

_ 

counter,  current 

! 

_ 1 

3. 1.2.5  INTERP 


Subroutine  INTERP  interpolates  between  tabulated  acceleration  values  for 
the  acceleration  at  time  TT.  A  trapazoidal  interpolation  scheme  is  used. 


CALLING  STATEMENT  ARGUMENT  LIST 

Mnemonic 

Symbol 
Ref.  (1) 

Type 

IREF 

1*4 

Acceleration  table  acceleration  values 

TT 

t 

R*8 

Integration  time  value 

YOUT 

R*8 

J 

Interpolated  acceleration  value 
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3. 1.2. 6  DERIV 


In  subroutine  DERIV  the  equations  of  motion  associated  with  each  mass 
element  are  formed.  The  terms  included  in  the  equations  are  based  on  the 
linear  and  nonlinear  stiffness  characteristics  and  the  damping  properties 
of  the  beam  elements.  Section  1.3. 5. 3,  Reference  (1);  the  external  aero¬ 
dynamic  forces  acting  on  the  vehicle,  Section  1.3. 5. 2,  Reference  (1);  the 
external  spring  forces  acting  on  the  vehicle.  Section  1.3. 5. 4,  Reference  (1); 
and  Inertia  properties  of  the  mass  elements.  The  equations  of  motion  are 
integrated  numerically.  Section  1.3.13,  Reference  (1),  to  determine  the 
incremental  displacement,  velocity,  and  acceleration  of  each  mass  element. 

External  spring  forces  and  deflections  are  supplied  to  DERIV  by  a  call 
to  subroutine  CFORCE.  Subroutine  ACCELT  is  called  only  when  the  input  option 
to  supply  time-history  tables  for  mass  accelerations  is  exercised.  All  data 
input  and  output  to  DERIV  is  accomplished  through  the  various  common  block 
regions  specified  in  the  subroutines. 

3. 1.2.7  CFORCE 

In  subroutine  CFORCE  the  computations  corresponding  to  the  theoretical 
analyses  discussed  in  Section  1.3. 5. 4,  Reference  (1),  are  carried  out.  The 
deflections  of  the  external  springs  which  are  in  contact  with  the  ground 
plane  are  determined.  Based  on  these  deflections,  the  force  acting  along 
the  axis  of  each  spring  is  determined.  The  component  of  this  force  acting 
perpendicular  to  the  ground  plane  is  combined  with  the  friction  coefficient 
supplied  in  the  input  data  to  give  the  drag  force  acting  on  the  spring. 

The  force  in  the  spring  and  the  drag  force  are  resolved  into  forces  and 
moments  acting  on  the  mass  element  to  which  the  spring  is  attached.  These 
forces  and  moments  are  returned  to  DERIV,  the  calling  subroutine,  through 
the  various  common  block  regions  included  in  the  subroutines. 

3. 1.2.8  INPUT 

Subroutine  INPUT  reads  all  input  card  data.  The  data  are  processed  and 
stored  in  appropriate  common  block  regions.  A  detailed  description  of  all 
input  data  required  by  program  KRASH  is  presented  in  Reference  (2). 
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3. 1.2.9  INPRNT 


Subroutine  INPRNT  processes  and  prints  data  supplied  by  subroutine  INPUT 
and  when  appropriate,  subroutine  GENMOD.  The  data  is  supplied  through  the 
various  common  block  regions  included  in  the  subroutine.  The  data  printed 
by  this  subroutine  are  discussed  in  detail  in  Reference  (2)  and  an  example 
listing  of  the  data  is  shown  in  Figure  3-3. 

3.1.2.10  ECHO 

Subroutine  ECHO  prints  a  complete  listing  of  the  input  data  in  card 
image  format.  An  example  of  this  output  is  shown  in  Figure  3-4. 

3.1.2.11  GENMOD 

Subroutine  GENMOD  generates  the  necessary  input  data  to  expand  a  half 
airplane  model  representation  into  a  full  airplane  model  representation. 

This  subroutine  is  called  only  when  the  RUNMOD  flag  is  set  equal  to  or  greater 
than  2.  The  various  common  block  regions  included  in  GENMOD  are  used  for  the 
input  and  output  of  data. 

3.1.2.12  RSOUT 

In  subroutine  RSOUT  all  state  vector  data  are  processed  for  subsequent 
output  through  calls  to  subroutine  DATOUT.  RSOUT  is  called  only  when  one  of 
the  check-point  times  of  the  restart  option  has  been  encountered.  Input 
to  the  subroutine  is  made  through  the  included  common  block  regions  and  the 
subroutine  call  argument:  list. 


CALLING  STATEMENT  ARGUMENT  LIST 

Symbol 

Mnemonic 

Ref.  (1) 

Type 

MODEL 

Name 

R*8 

Check  point  run  name  identifier 

CASE 

Case  No. 

1*4 

Check  point  case  number  identifier 

MSEC 

1*4 

Check  point  time  in  milliseconds 
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Figure  3-3.  Example  INPRNT  Output  (Sheet 


VEHICLE  ROTATIONAL  VflOCITUS  IN  VEHICLE  AXIS  I  RAO/SEC  I 
EULER  ANGLES  OF  VEHICLE  RELATIVE  TO  GROUtO  I  RADIANS  I 


Figure  3-3.  Example  INPRNT  Output  (Sheet 


MtlGHTS  HASS  COORDINATES  F.S. >B.L. >H.L.  HASS  HOMCNTS  OF  INERTIA  4 IB-IH-JEC»»E I 


•  0000000000000000000000-^00000000 
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3.1.2.13  RSIN 


In  subroutine  RSIN  state  vector  data  are  retrieved  through  calls  to 
subroutine  DATIN.  RSIN  is  called  only  when  a  restart  run  is  made.  Input 
of  data  to  the  subroutine  is  made  through  the  calling  statement  argument 
list  while  output  is  through  the  included  common  block  region. 


CALLING  STATEMENT  ARGUMENT  LIST 

Symbol 

Mnemonic 

Ref.  (1) 

Type 

IMODEL 

Name 

mm 

Restart  run  name  identifier 

ICASE 

Case  No. 

Restart  run  name  identifier 

IMSEC 

19 

Restart  time  in  milliseconds 
_ * 

3.1.2.14  DOAIJ 

Beam  fundamental  uncoupled  frequencies,  beam  constants,  and  beam  Euler 
angle  transformations  are  computed  in  subroutine  DOAIJ.  The  frequencies  are 
computed  using  equations  l-55a  and  l-55b  of  Reference  (1)  while  the  damping 
constants  are  determined  using  equation  1-54  of  Reference  (1).  All  input  and 
output  to  the  subroutine  is  accomplished  through  the  various  common  block 
regions  included  in  the  subroutine. 

3.1.2.15  RC 

Within  subroutine  RC  the  initial  conditions  for  all  model  elements  are 
computed  for  a  restart  run.  In  addition,  allowable  force  and  deflection 
are  determined  and  printed.  Subroutine  RC  is  called  only  when  a  restart  case 
is  run.  Input  and  output  to  the  subroutine  is  accomplished  through  the  vari¬ 
ous  common  block  regions  included  in  the  subroutine. 


3.1.2.16  PRINT 


Subroutine  PRINT  is  called  at  the  beginning  of  each  run  and  at  each 
print-time  interval  specified  in  the  input  data.  The  subroutine  provides 
for  the  printing  of  dynamic  response  information  defining  the  motion  of 
massless  nodes  and  masses,  the  forces  and  deflections  associated  with  all 
beam  and  spring  elements,  and  a  complete  energy  distribution  summary  for 
the  model.  Examples  of  the  data  printed  are  shown  in  Figure  3-5. 

When  summary  plots  are  requested,  the  necessary  mass  and  element  data 
are  processed  and  written  on  the  temporary  data  storage  units  specified  in 
Section  2.2.  Storage  of  this  data  takes  place  at  each  print-time  interval. 

Position  plots  are  also  generated  within  the  PRINT  subroutine  by  a  call 
to  subroutine  PAPLOT.  This  call  is  made  at  the  plot  time  interval  specified 
in  the  input  data. 

All  input  and  output  to  subroutine  PRINT  is  accomplished  through  the 
various  common  block  regions  included  in  the  subroutine. 

3.1.2.17  IC 

Within  subroutine  IC  the  initial  condition  data  for  all  elements  of 
the  model  are  determined.  Subroutine  IC  is  called  at  the  beginning  of  each 
run  except  a  restart  run,  in  which  case  subroutine  RC  is  called.  In  addi¬ 
tion  to  the  initial  condition  data,  subroutine  IC  determines  allowable  beam 
force  and  deflection  data.  These  data  are  based  on  the  theoretical  discus¬ 
sion  presented  in  Section  1.3.17  of  Reference  (1). 

An  example  of  the  data  printed  from  subroutine  IC  is  shown  in  Fig¬ 
ure  3-6.  Input  and  output  to  this  subroutine  is  accomplished  through  the 
various  common  block  regions  included  in  the  subroutine. 
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Figure  3-5.  Example  PRINT  Output  (Sheet 


SIAM  F04C  FS 


©  ©  ©  —  o  c  oc»  oooo 

OC©—  COCC©  coo*- 


2  «©»  ©—©—©©—*©*  °. 

^  u,  r  «  —  «  or-©©©'©©** 

-  X*.—  —  I  «  l  I  c  i  I  l  l 


S£S©£ooo©  5  3  o  o 

“♦cooccoioc  cccc 

**©•*•©—  —  ©•*<  —  ©  —  ©  •- 


►  -  —  ©Or-  *  ^  ©  ' 

►  •  «>  ©  0  ©  ©  —  ►  • 


SoSoScooo  O  C  ©  © 

—  ©ooooooc©  e  ©  ©  c 

•W  0  fN.  —  ©  »“  —  «  0  —  ©  *■  ©  »- 

X  «©WN-W(rtr  ►-  0  O  © 

►  000©»C©e«  «.  —  r- 

«N4wr«i»iC  ©  ►  ©  • 


- -  _  o  *  r>  O  O  C  —  — 

cccoo©©©©  o  ©  c  © 

■>c:iic.c{!-;  c  c  j  c 

►  «*©*v. 

I  ^  f  —  r.0*  —  0 

♦  o 

**•©»>■»'.*»©>•*  0  «  *•  0 


»  -  »NNC  J  O.  —  O  ©  ©.  — 

ccocccc&c  O  ©  —  © 
—  c.  c  cccccu  cocc 

►  —  C  —COO  C  0  — 

IC’-r«ir»t* - o  «r  *  r 

»»e  "  f  c  ©  —  *  «r  —  0  c 

—  OOIN'  4  C  t  Q 


■fNNNO«"KP> 

-  6  ©  c  —  -  c  ©c< 
1*1111  I 
-CCCCCCCCf 

‘*•**©*■0*  IN  *V  , 

■  ►©  «  *  r>.  ©  *«•  0  P 
'  *  N  r  »  »  c  ^  »  I 
:»0«0N<N*< 


See-?©©  —  —  ©  cccc 
COCCCOCCC  cccc 

•*  «  ©  —  «M  —  ©  0  J  ©  *-  0  ©  «S. 

©00©©-*0«0  0—00 
~0«*  —  r»*rr*©®«  —  0 


oScooScec 

-»ce©eoec©c 

—  ©■'►**—©-©0«f'0 

X  o  *•  &  *  *  c  ••■  ©  c 


&o©o§»§>Soe 


L  r  «  ©  »  0  —  -  ©►• 

r*  ©*  r»  »  ©  0  o  e  0 

—  O  0  ar  0  «r  —  ^  « r 


—  —  —  —  —  0"'f.0 

ccocceocc 
Tiii  •  T 
-»cccccccdc 
►  *  *  ^  ^  <  »  r  n  fr 

r  ©  *■  pw  •&  *.  0  p-  »»  *•■ 

*s.©  ©  0  —  r.  f-  0  © 

0  IT  P-  0  —  ©■  P\.  P- 


-cccccccc 

X»M^«l^-*-» 
—  r  *■*  »r  r.  —  ►  >. 
P-  ©  0,  0  —  P-  * 


o  —  —  ©eo  —  —  £ 

I  I  I  I  I  I  I  I 

CCCCOC  DCC 

:  —  ©  p-  P*  *r  0  p.  c  0 
©P>Pw,P-«ew*s0^ 
P-0CO®©— «© 


CCOOCCOCC 
■wc-c  —  Or*-.*  — 
u.  ©  p»  p">  0  0  ©  ©>  *  0. 

»  NCK  C 

NN<*-OOI'*r 


c  c  c  c  c  e  < 

—  —  ©  0  ^  ©  0  0  * 

—  ^0p’-rv.«if  —  • 

*  k  <  n*.  ►  r  ►  c 

-">  0  p-  0  ^  O  • 


p-  m  C  C  C  C  C  C  C 

c  —  ►<•*«©©»< 

0  —  C  0  —  —  rvf 

2  Sc  S££3; 


PS  P  rvlf  0  0  iir  N*  < 

77T<:cr<r<T<f??<T<f<f 

ccccecccccc  e© 

—  «r  «•.«.  —  p.  —  p  A  —  *'■'©»*• 
«*  —  e0p*Pve‘cC*-p'0  — 
a  «  c  r*  p  0  0  ©  0  p  O  —  «■  «f 
►  P.  ©  ©  p-  0  C  p.  0  0  0  0  0 


p  J  -  j  •«  w  0  P*  0  p*  *  V 

cceeoocccooc© 

—  ©  "  *-  0  *^  —  <■.  00  •*  *  •• 

u,  p-  *-  «  00&P-—  0'0*'<-'‘© 

z  0  «r  r  e  00^  —  oc  «  ;  c 

—  ,r  —  <*0f^0*»>-c«rp‘C  © 


— .  p«le>(rir00»'.Pir©p»i»p 

— - ecc  —  —  c  ©occ 


CCCCCCCCCOCCC 
»»'fC^^-irN»  C  ©  « 
.*■  if  t  *■  »  r—  ©r^f'C^ 
00  fs.a.  —  p*.  0  c  —  ©  0  ©  0 


C  ©  —  ©  —  *M  e*  • 
- - S  O  C  -  • 

I  I  I  I  ^  I  I 

>.  t  Sr. 

©  ©  ©  p*  p>  0,  —  • 

—  f  p-  r*  "  **  ©  • 


C  ©  C  C  C  C  ccc 

►  O  —  ©»r  —  0  &  p.  r* 

U  4  —  P-  *N<  tr  ©«c*r' 

—  0*00rp-O  — 
©  p*  ©  0  O  ©  k  P’  0' 


E.;r?rcff:prif? 

-  .f  "  »^IN.--~*S.~©<rf 

©  -■  .  kf  -V.  fNv  —  ».  Cr  — 

■»  j-  -r  —  c  f'  r  rtf. 


—  C  ©  —  C  C  « 

©  C  O  O  ©  C  < 

?  C  C  c  O  C  £ 

k  ©  <n>  r  ^  *  *-  c 

U.f^©  »«■••<( 

•r  *  r  o  —  m 

»*l".  Ol*  Nf 


—  ©  ©  —  ©  C  —  *M 

©C  CONOCO 

ic  occocc 

k  <  n  r  n  n  e  i*  c 

U.  —  ^  —  —  00  —  * 

e*  ir  c  if  c  c  •■ 

fw-COiNf'f.  * 


N  N  c  f*  <  I 


ec-cceccccocc 

(0—  ©>«^0«rC  /  f'fs.t 
—  C-' 

c  —  r-f*  — O"  »-  c©*mC*v- 


xoccooocc  © 

»  xeeojeococ 

* 

M.  *»NF  mf  ■#  IT4  *  ►  i 


rocoooocc© 


X  ©OOOOSOOC 


XCOC©OCC©C—  -CO 


r  croorrc-'-'Cc  rc  © 


^^»rr*lr*<f‘fc|f'^l, 


•  *v  —  iv  — 


Figure  3-5.  Example  PRINT  Output  (Sheet  2  of  4) 


I  IN  Ota  jdllHiil  «|  A *41  V  »VM 


Figure  3-5.  Example  PRINT  Output  (Sheet  3  of  4) 


y  oeoooooc 
k  oooocooc 


»  © —00  —  00 


Uk.  «  _ 

“  *  ooooorj 


£  0000000c 

00—00—00 


D  ©  O  o  ©  O  ©  ( 

*  O  O  O  O  o  O  < 


Z  Z  o  O  «  o  ©  , 


OO^CC*  00 
©©Voodoo 


£  >  c  v 
—  c  r  * 


*  ©  O  o  O  c  O  < 


^C^OC-^CC 


£  ©©©OCCO—Cfw 

a  ccc-dococ  —  c 


30  —  00  —  cc 


*>  ©ooc-eocooc 
•  Sr 

J  5  ^  e.  *,  'J  r,  N.  " 


£  ®®©©ooooco 
*  CCOCCOOCNC 


CO«oo«CO 

OOiroorcc 


oo*co«rcc 

60-00-00 


C“*“fV—  CN. 

*  ococcooooc 

J1  *  OCOOOOOOOO 

X  z 

Bt«  ■»  r*~*'  *■  +  *■■***.*. 


oocococs  c 
©OOCoOCO  £ 


©  ©  »  oo*^oc  v  M 1 
oo-do-eo  -  J  * 

ceoocoec  ^ 


—  —  •»*«*-  -  *-  r* 

- E  S- 


t  >-■  ecc.c  =  cceec 

-SffEjsfls 


C  5  *  ‘  *  «  «  =  »  c 

6  r  -  -  N  <  N  <  f  « 


—  Z  S  P.  _ 

-  :  ► 


►  er'  —  r 


■  5  .*..*,*  0  *  *•  * 

*-  *  f  *■’  —  o*  -  J  ^ 


Figure  3-5.  Example  PRINT  Output  (Sheet  4  of  4) 


Figure  3-6.  Example  IC  Output  (Sheet  1  of  2) 


•fan  r*>*i  rr  ti>ns 


Figure  3-6.  Example  IC  Output  (Sheet  2  of  2) 


3.1.2.18  MATMUL 


Subroutine  MATMUL  carries  out  the  matrix  multiplication  A*B  and  returns 

the  results  in  C. 


CALLING  STATEMENT  ARGUMENT  LIST 


Mnemonic 

Symbol 
Ref.  (1) 

Type 

A 

R*8 

Input  matrix 

B 

R*8 

Input  matrix 

C 

R*8 

Resulting  Output  matrix 

3.1.2.19  PREPLT 

Subroutine  PREPLT  is  called  at  the  end  of  each  run.  The  subroutine 
retrieves  and  reorders  data  sets  previously  stored  on  the  logical  units 
specified  in  Section  2.2.  After  processing,  the  data  are  transferred  to 
subroutine  PLOTT  for  subsequent  output  as  summary  plots.  When  the  processing 
of  data  cannot  be  accomplished  within  the  working  space  allocated  to  the  sub¬ 
routine,  logical  unit  12  is  used  for  the  temporary  storage  of  the  processed  data 
The  maximum  number  of  elements  that  can  be  accommodated  in  the  allocated  working 
space  is  tabulated  in  Table  3-2.  Input  of  data  to  the  subroutine  is  through  the 
calling  statement  argument  list  as  well  as  through  the  included  common  block 
regions. 


CALLING  STATEMENT  ARGUMENT  LIST 

Mnemonic 

Symbol 
Ref.  (1) 

Type 

DELTAT 

Integration  time  interval 

NM 

1 

Number  of  mass  elements 

NNP 

1*4 

Number  of  massless  node  elements 

1*4 

Number  of  beam  elements 

1*4 

Number  of  external  springs 

NDRI 

1*4 

Number  of  DRI  elements 

MSECIN 

1*4 

Check  point  save  times 

TABLE  3-2.  PLOT  ELEMENT  WORKING  SPACE  LIMITS 


Element 

Maximum  Number  of  Elements 

Accommodated  in  Working  Space 

MASS 

31 

MASSLESS  NODE 

41 

BEAM  (Forces) 

51 

BEAM  (Deflections) 

51 

BEAM  (Stresses) 

51 

EXTERNAL  SPRINGS 

ul 

BEAM  (Energies) 

51 

DRI 

101 

3.1.2.20  PLOTT 


Subroutine  PLOTT  is  called  by  subroutine  PREPLOT  and  provides  printer 
plotting  capability  for  summary  time  history  plots  of  requested  model  elements. 
This  subroutine  is  a  modified  version  of  one  published  by  Dr.  David  H.  Laananen, 
Reference  (4).  The  subroutine  accepts  data  through  the  calling  statement 
argument  list  and  outputs  the  data  in  a  special  plot  format.  An  example  of  the 
plot  format  is  shown  in  Figure  3-7. 


CALLING  STATEMENT  ARGUMENT  LIST 

Mnemonic 

Symbol 

Ref.  (1) 

Type 

M 

1*4 

Number  dependent  variables 

NP 

1*4 

Dimension  of  X,  Yl,  Y2 ,  Y3  arrays 

X 

R*4 

Values  for  the  independent  variable 

Y1 

R*4 

) 

Y2 

R*4 

\  Values  for  the  dependent  variables 

Y3 

R*4 

1 
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Figure  3-7.  Example  PLOTT  Output 


3.1.2.21  PAPLOT 


Subroutine  PAPLOT  is  called  by  subroutin  PRINT  at  the  beginning  of  each 
ran  and  at  each  position  plot  time.  See  Section  2.2.7,  Reference  (2).  The 
subroutine  provides  program  KRASH  with  on-line  plotting  capability  for  plotting 
the  location  of  mass  elements  in  two-dimendional  space.  An  example  plot  is 
shown  in  Figure  3-8. 

Calls  are  made  from  PAPLOT  to  subroutines  FSHELL,  SHELLX,  and  SHELLM 
where  data  to  be  plotted  is  sorted.  Based  on  the  returned  sorted  data,  format 
statements  are  generated  which  produce  the  desired  position  plots.  Input  data 
is  supplied  to  PAPLOT  through  the  calling  statement  argument  list  and  the 
included  common  block  regions. 


CALLING  STATEMENT  ARGUMENT  LIST 

Mnemonic 

Symbol 
Ref.  (1) 

Type 

V 

R*8 

Quantity  to  be  plotted  along  horizontal 
axis 

W 

R*8 

Quantity  to  be  plotted  along  vertical 
axis 

XSCALE 

R*4 

Horizontal  axis  scale  factor 

Y SCALE 

R*4 

Vertical  axis  scale  factor 

NPTS 

1*2 

Number  of  point  to  be  included  in  plot 

ISCALE 

1*2 

Scale  option  control  flag 
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Figure  3-8.  Example  PAPLOT  Output 


3.1.2.22  FSHELL 


Subroutine  FSHELL  is  called  by  subroutine  PAPLOT.  In  this  subroutine 
elements  located  in  ARRAY  are  sorted  and  relocated  in  ARRAY  in  ascending  order 
(algebraically  least  first).  The  array  KEY  identifies  and  returns  the  original 
location  of  the  corresponding  element  in  ARRAY. 

All  data  are  passed  through  the  calling  statement  argument  list. 


CALLING  STATEMENT  ARGUMENT  LIST 


Mnemonic 

Symbol 
Ref.  (1) 

Type 

ARRAY 

R*4 

Array  containing  elements  to  be  sorted 

KEY 

1*4 

Array  containing  integer  flags  which 
identify  the  location  of  corresponding 
ARRAY  elements  prior  to  sorting 

N 

1*4 

Number  of  elements  in  ARRAY 

3.1.2.23  SHELLX 


Subroutine  SHELLX  is  called  by  subroutine  PAPLOT.  In  this  subroutine 
elements  of  the  array  DARRAY  are  reordered  in  accordance  with  the  locations 
specified  by  the  elements  in  KEY. 


All  data  are  passed  through  the  calling  statement  argument  list . 


CALLING  STATEMENT  ARGUMENT  LIST 

i 

Mnemonic 

Symbol 
Ref.  (1) 

Type 

DARRAY 

R*4 

Array  containing  elements  to  be  sorted 

KEY 

1*4 

Array  containing  integer  flags  which 
identify  order  in  which  elements  of 

DARRAY  are  to  be  sorted. 

N 

1*4 

Number  of  elements  in  DARRAY 

3.1.2.24  SHELLM 


Subroutine  SHELLM  is  called  by  subroutine  PAPLOT.  In  this  subroutine 
elements  of  the  array  M  are  reordered  in  accordance  with  the  locations  specified 
by  the  elements  in  KEY. 

All  data  are  passed  through  the  calling  statement  argument  list. 


CALLING  STATEMENT  ARGUMENT  LIST 

Mnemonic 

Symbol 
Ref.  (1) 

Type 

KEY 

1*4 

Array  containing  integer  flags  which 
identify  order  in  which  elements  of  M 
are  to  be  sorted. 

M 

1*4 

Array  containing  elements  to  be  sorted 

N 

1*4 

Number  of  elements  in  M 

3.2  COMMON  BLOCK  REGIONS 

3.2.1  Common  Block/Subroutine  Cross  Reference 

Program  KRASH  is  designed  so  that  data  storage  and  transfer  is  accomplished 
using  the  many  common  block  regions  defined  within  the  program.  A  cross 
reference  of  the  common  block  names  and  using  subroutines  is  given  in  Table  3-3. 
Included  in  the  cross  reference  summary  are  the  core  sizes  required  by  the 
IBM  360  computer  system  for  each  of  the  regions. 

3.2.2  Symbol  Definitions 

All  of  the  common  block  regions  allocated  in  program  KRASH  are  summarized 
in  Table  3-4.  Included  in  the  summary  are  the  symbols  used  to  identify  the 
variables 'which  make-up  each  block  along  with  the  symbol  type,  array  dimension, 
and  a  brief  definition. 
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Subroutine  FSHELL  is  called  by  subroutine  PAPLOT.  In  this  subroutine 
elements  located  in  ARRAY  are  sorted  and  relocated  in  ARRAY  in  ascending  order 
(algebraically  least  first).  The  array  KEY  identifies  and  returns  the  original 
location  of  the  corresponding  element  in  ARRAY. 

All  data  are  passed  through  the  calling  statement  argument  list. 


CALLING  STATEMENT  ARGUMENT  LIST 


Mnemonic 


Symbol 
Ref.  (1) 


Type 


ARRAY 

KEY 


R*4 

1*4 

1*4 


Array  containing  elements  to  be  sorted 

Array  containing  integer  flags  which 
identify  the  location  of  corresponding 
ARRAY  elements  prior  to  sorting 

Number  of  elements  in  ARRAY 


3.1.2.23  SHELLX 

Subroutine  SHELLX  is  called  by  subroutine  PAPLOT.  In  this  subroutine 
elements  of  the  array  DARRAY  are  reordered  in  accordance  with  the  locations 
specified  by  the  elements  in  KEY. 

All  data  are  passed  through  the  calling  statement  argument  list. 


CALLING  STATEMENT  ARGUMENT  LIST 

Mnemonic 

Symbol 
Ref.  (1) 

Type 

DARRAY 

R*4 

Array  containing  elements  to  be  sorted 

KEY 

1*4 

Array  containing  integer  flags  which 
identify  order  in  which  elements  of 

DARRAY  are  to  be  sorted. 

N 

1*4 

Number  of  elements  in  DARRAY 
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TABLE  3-3.  COMMON  BLOCK/SUBROUTINE  CROSS  REFERENCE  SUMMARY 


Si  ze  in 

By  tes 
Fortran  H 
Compiler 
(OPT  =  2) 

Subrout ines 

ACCELT 

CDATE 

CFORCE 

DAT  IN 

DATOUT 

DERIV 

DOAIJ 

ECHO 

EULER 

FSHELL 

GENMOD 

o 

►-< 

INPRNT 

INPUT 

INTERP 

JOBAC 

MAIN 

MATMUL 
MAT V EC 
PAPLOT 
PLOTT 

PREPLT 

PRINT 

RC 

RSIN 

RSOUT 

BLANKCOM 

9980 

• 

• 

• 

• 

CFIC 

58 

• 

• 

• 

CFIR 

640 

• 

• 

CFPR 

B4D 

• 

• 

• 

• 

COMALL 

1C088 

• 

•  • 

• 

• 

• 

• 

• 

• 

•  • 

• 

COMNEW 

18 

• 

• 

• 

• 

COMMI2 

8D4 

• 

• 

• 

COMMI4 

268 

• 

• 

• 

• 

• 

COMMR4 

2FDO 

• 

• 

• 

COMMR8 

4B0 

• 

• 

• 

DEIC 

260 

• 

• 

• 

• 

DEIN 

1022C 

• 

• 

• 

• 

• 

• 

» 

DEIR 

1 2CO 

• 

• 

DEINPR 

48B0 

• 

• 

• 

• 

• 

• 

• 

DEMA 

8 

• 

• 

DEPR 

67B0 

• 

• 

• 

• 

DINICP 

2044 

• 

• 

• 

• 

• 

• 

•  • 

• 

DOIN 

4B0 

• 

• 

• 

• 

ENERGY 

16B  8 

• 

• 

• 

• 

I  BALL 

50 

• 

• 

INAC 

1 5E4 

• 

• 

• 

INCF 

DC  8 

• 

• 

• 

• 

INCFIC 

8 

• 

• 

• 

• 

• 

• 

INDEAC 

4 

• 

• 

• 

• 

• 

INIC 

CD8 

• 

• 

# 

• 

• 

• 

INIDCP 

280 

• 

•  • 

• 

• 

• 

• 

• 

• 

ININPR 

20 

• 

• 

• 

INOUT 

Li _ 

• 

• 

• 
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RSIN 

RSOUT 

SHELLM 

SHELLX 


TABLE  3-3.  COMMON  BLOCK/ SUBROUTINE  CROSS  REFERENCE  SUMMARY  (Continued) 


Common' 

Block 

Regions 

Fortran  H 
Compiler 
(OPT  -  2) 

Subroutines 

ACCELT 

CDATE 

CFORCE 

DATIN 

DATOUT 

DERIV 

DOAIJ 

ECHO 

EULER 

FSHELL 

GENMOD 

u 

►-H 

INPRNT 

INPUT 

INTERP 

JOBAC 

MAIN 

MATMUL 

MAT  V  EC 

PAPLOT 

PLOTT 

PREPLT 

PRINT 

o 

oc 

RSIN 

RSOUT 

5  3 

J  -2 

U-2  UJ 

X  X 
in  c n 

INPR 

8 

• 

• 

• 

• 

• 

•  • 

• 

• 

• 

IPIC 

4 

• 

• 

• 

• 

IRDE 

8 

• 

• 

MACF 

3E84 

• 

• 

• 

• 

• 

MACFIN 

C 

• 

• 

• 

• 

•  • 

• 

MADE 

1134 

• 

• 

• 

• 

MAPR 

12CC 

• 

• 

• 

• 

MAPR1 

976 

• 

• 

MAX 

2B98 

• 

# 

• 

• 

MCFIII 

8 

• 

• 

• 

• 

•  • 

• 

• 

NP00I2 

320 

• 

•  « 

• 

• 

• 

•  • 

• 

• 

NP00I4 

4 

• 

•  • 

• 

• 

• 

•  • 

• 

• 

• 

• 

NP00R8 

4B0 

• 

• 

• 

• 

NP01I2 

F0 

• 

• 

• 

• 

• 

• 

• 

NP01R8 

1C20 

•  • 

• 

• 

• 

• 

• 

NP0212 

780 

• 

• 

• 

NP02R8 

4B0 

• 

• 

• 

• 

• 

• 

• 

OLEO 

82C 

• 

• 

• 

• 

• 

• 

• 

OTPLT 

F72 

•  • 

• 

• 

PPLTS 

4CA 

• 

• 

• 

PRMA 

2490 

• 

• 

• 

• 

• 

RESTRT 

30 

• 

•  • 

STUFF 

3A0 

• 

• 

UB 

70C 

• 

• 

• 

• 

VARINT 

40 

_ _ 

• 

• 

•  • 

• 
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COMMON  BLOCK  SUMMARY 
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TABLE  3-4.  COMMON  BLOCK  SUMMARY  (Continued) 


Common 

Name 

Mnemonic 

Type 

Array 

Size 

BLANKCOM 
(Cont' d) 

WOLD 

R*8 

80 

XOLD 

R*8 

80 

YOLD 

R*8 

80 

ZOLD 

R*8 

80 

PINO 

R*8 

80 

QINO 

R*8 

80 

RINO 

R*8 

80 

PHIOLD 

R*8 

80 

THEOLD 

R*8 

80 

PSIOLD 

R*8 

80 

XACFD 

R*8 

80 

YACFD 

R*8 

80 

ZACFD 

R*8 

80 

XAFOLD 

R*8 

80 

YAFOLD 

R*8 

80 

ZAFOLD 

R*8 

80 

XNPFD 

R*8 

SO 

R*8 

50 

ZNPFD 

R*8 

50 

XANPFO 

R*8 

50 

YANPFO 

R*8 

50 

ZANPFO 

R*8 

50 

TKR 

R*8 

TPEN 

R*8 

80 

DTHALF 

R*8 

Previous 
Value  for 


Previous 
Value  for 


Time 

Derivative 

of 


Previous 
Value  for 


Time 

Derivative 

of 


Previous 
Value  for 

Times  for  beam  yielding  and 
rupture  summary  table. 

Penetration  times. 

Half  of  time  interval. 
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TABLE  3-4.  COMMON  BLOCK  SUMMARY  (Continued) 
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TABLE  3-4.  COMMON  BLOCK  SUMMARY  (Continued) 


Common 

Array 

Name 

Mnemonic 

Type 

Size 

COMALL 

VEE 

R*8 

900 

Beam  relative  deflection 

(Cont' d) 

vector . 

■ 

R*6 

80 

Mass  weight. 

R*8 

80 

)  Euler  angles  defining 

R*8 

80 

>  orientation  of 

mass  I 

PSI 

R*8 

80 

^  relative  to  ground. 

PDOT 

R*8 

80 

\  i 

1  P 

QDOT 

R*8 

80 

Q 

R*8 

80 

I 

R 

R*8 

80 

1 

u 

R*8 

80 

I 

V 

WDOT 

R*8 

80 

1  Time 

\  derivative 

/  w 

XDOT 

R*8 

80 

/  of 

\  x 

YDOT 

R*8 

80 

i 

1  Y 

ZDOT 

R*8 

80 

\ 

f  7 

PHIDOT 

R*8 

80 

I 

PHI 

THEDOT 

R*8 

80 

THETA 

PSIDOT 

R*8 

80 

/ 

\  PSI 

TIME 

R*8 

Time . 

DFLTAT 

R*8 

Integration  time 

interval . 

XACC 

R*8 

80 

YACC 

R*8 

o 

00 

{  Acceleration  in 

(  body  axis. 

ZACC 

R*8 

80 

) 

AIDOT 

R*8 

9 

PHIIJ 

R*8 

150 

l  Euler  angles  defining 

THEIJ 

R*8 

150 

(  beam  axis  relative  to 

i  ground.  Used  to  calculate 

PSIIJ 

R*8 

150 

;  Au. 

SUMDF 

R*8 

6,150 

Summation  of  incremental 
strain  forces  over  time. 
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TABLE  3-4.  COMMON  BLOCK  SUMMARY  (Continued) 


Common 

Array 

Name 

Mnemonic 

Type 

Size 

COMALL 

TITLE 

1*4 

40 

Run  title. 

(Cont' d) 

XLBAR 

R*8 

40 

External  spring  length. 

FSPBAR 

R*8 

40 

External  spring  axial 

load 

when  unloading  occurred 

• 

VEEDOT 

R*8 

3.  3 

Spring  contact  point 

velocities  relative  to 

ground. 

DX 

R*8 

81 

)  (  x 

DY 

R*8 

81 

f  change,  over  ) 

/  time  of  S 

DZ 

R*8 

81 

)  (  z 

DP  IN 

R*8 

81 

)  change  (  vP ^  )  in< 

J  -  us:  I™ 

DQIN 

R*e 

81 

:remental 

tation 

DRIN 

R*8 

81 

SEIJ 

R*8 

150 

Strain  energy  for  beam 

IJ. 

DEIJ 

R*8 

150 

Damping  energy  for  beam 

IJ. 

CEIK 

R*8 

40 

Crushing  energy 

for  ex- 

CEIKF 

R*8 

40 

Friction  energy 

ternal 

SBAR 

R*8 

40 

Unloading  deflection 

spring 

IK 

KUN 

R*8 

40 

Unloading  stiffness 

MAXNM 

1*4 

Maximum  number  of  masses. 

MAXIGS 

1*4 

Maximum  number  of  beams 

• 

MAXTBL 

1*4 

Maximum  number  of  KR 

tables . 

NM 

1*4 

Number  of  mass  elements 

• 

NB 

1*4 

Number  of  beam  elements 

• 

I 

1*4 

Mass  number. 

J 

1*4 

Mass  number. 

IG 

1*2 

150 

Mass  at  ic^  ^  end 

of  beam 

JG 

1*2 

150 

Mass  at  j1'1  j  IJ 

beam 

spring 

N1 

1*2 

900 

Unloading  -  reloading 

NN 

1*2 

40 

Flag  for 

IJPR 

H 

150 

DRI  element  flag 
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TABLE  3-4.  COMMON  BLOCK  SUMMARY  (Continued) 


Common 

Name 

Mnemonic 

Type 

Array 

Size 

COMM  12 

IQ 

1*2 

180 

Mass  number  1  for  KR  table. 

JQ 

1*2 

180 

Mass  number  J  of  KR  table. 

LQ 

1*2 

180 

Direction  identification  for  KR  table. 

NPQ 

1*2 

180 

No.  of  points  in  KR  table. 

INBUF 

1*2 

80 

Mass  numbers  for  nonstandard  HE's. 

NKMVEC 

1*2 

150 

Beams  with  directly  input  K  matrices. 

ISAVE 

1*2 

180 

Beams  with  KR  tables. 

COMMI4 

ND 

1*4 

No.  of  nonstandard  damping  values  for 
beams. 

NVBMN 

1*4 

No.  of  beams  having  nonstandard  maximum 
negative  deflection  cut-off. 

NFBMN 

1*4 

No.  of  beams  having  nonstandard  maximum 
negative  load  cut-off. 

NHI 

1*4 

No  of  nonstandard  He's  or  non  zero  IxY> 

iyz  >  izx- 

NKM 

1*4 

No.  of  user  supplied  stiffness  matrices. 

NLB 

1*4 

No  or  nonlinear  beam  KR  curves. 

NPH 

1*4 

No.  of  masses  having  nonstandard  Euler 
angles . 

NMTL 

1*4 

No.  of  nonstandard  materials. 

NPTS 

1*4 

50 

No.  of  points  in  accel.  history  table. 

NVBM 

1*4 

No.  of  beams  having  nonstandard  maximum 
positive  deflection  cut-off. 

IJPRT 

1*4 

14 

DR1  element  identification  numbers. 

IPHDP 

1*4 

80 

Mass  number  for  nonzero  Euler  angles. 

NFBM 

1*4 

No.  of  beams  having  nonstandard  maximum 
positive  load  cut-off. 

COMMR4 

KR 

R*4 

2  700 

Stiffness  reduction  factor. 

LDP 

R*4 

180 

Linear  deflection  point. 

LDP  1 

R*4 

180 

Linear  deflection  point  for  reloading. 

COMMR8 

G 

R*8 

150 

Beam  material  shear  modulus. 

COMNEW 

DAMPC 

R*8 

Default  critical  damping  coefficient  for 
all  beams. 

RUNMOD 

R*8 

Run  mode  control  flag. 

PUNMOD 

R*8 

Punch  mode  control  flag 
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TABLE  3-4.  COMMON  BLOCK  SUMMARY  (Continued) 


Common 

Array 

-  -  , 

Name 

Mnemonic 

Type 

Size 

DEIC 

WTOT 

R*8 

Total  vehicle 

weight . 

CLTEST 

R*4 

150 

Beam  elements 

which  lie  in  the  airplane 

plane  of  symmetry. 

DEIN 

XNBAR 

R*8 

> 

XPBAR 

R*8 

YNBAR 

R*8 

Coordinates  defining  volume  for  MASS 
►  penetration  calculations. 

YPBAR 

R*8 

ZNBAR 

R*8 

ZPBAR 

R*8 

VOLENZ 

R*8 

5,  3 

Volume  change 
lengths . 

reference  volumes  time  zero 

HEX 

R*8 

80 

■N 

HEY 

R*8 

80 

-  Angular  momenta  of  rotating  masses. 

HEZ 

R*8 

80 

- 

ALIFT 

R*8 

80 

Aerodynamic  lift  constant  for  mass  I. 

VMAX 

R*8 

900 

Maximum  allowable  beam  positive 

deflections  - 

rupture. 

VMAXN 

R*8 

900 

Maximum  allowable  beam  negative 

* 

deflections  - 

rupture. 

FMAXN 

R*8 

900 

Maximum  allowable  beam  negative  force  - 

rupture . 

XKS 

R*4 

2700 

Slopes  ' 

XKI 

R*4 

2700 

Intercepts 

>-  for  KR  table  calculations. 

XKR 

R*4 

2700 

Deflections 

NLSFLG 

1*2 

900 

Nonlinear  beam  table  member. 

i 

CHUG 

1*2 

180 

5 

tarting  point  in  KR  table  for  beam  IJ. 

MVP 

1*4 

Mass  volume  penetration  identification. 

FMAX 

R*8 

900 

Maximum  allowable  beam  positive  force  - 

_ 

_ 

rupture . 
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TABLE  3-4.  COMMON  BLOCK  SUMMARY  (Continued) 


Common 

Array 

Name 

Mnemonic 

Type 

Size 

DEINPR 

AA 

R*8 

150 

Beam  cross  sectional  area. 

E 

R*8 

150 

Beam  material  elastic  mod. 

YY 

R*8 

150 

Area  moment  of  inertia  Iy. 

ZZ 

R*8 

150 

Area  moment  of  inertia  1%. 

XIQ 

R*8 

150 

Torsional  stress  parameter. 

XLB 

R*8 

150 

Beam  length. 

Z1 

R*8 

150 

Distance  to  extreme  fiber. 

Z2 

R*8 

150 

Distance  to  extreme  fiber. 

MC 

1*4 

150 

Material  code. 

XJ 

R*8 

150 

Beam  polar  area  moment  of  inertia. 

SF26 

R*8 

150 

Shape  factor  for  moment  about  Z  azis  at 
beam  ith  end. 

SF35 

R*8 

150 

Shape  factor  for  moment  about  Y  axis  at 
beam  i^h  end. 

SF26J 

R*8 

150 

Shape  factor  for  moment  about  Z  axis  at 
beam  jth  end. 

SF35J 

R*8 

150 

Shape  factor  for  moment  about  Y  axis  at 
beam  j ^  end. 

PY 

1*4 

150 

L  Designates  beam  pinned  at  ich  end. 

PZ 

1*4 

150 

J 

PYJ 

1*4 

150 

L  Designates  beam  pinned  at  j1-*1  end. 

PZJ 

1*4 

150 

/ 

NSC 

1*4 

Flag  to  call  for  stress  calculation. 

NPIN 

1*4 

150 

No.  of  pin-ended  beams. 

DEIR 

PLM26 

R*8 

150 

Plastic  moment  about  Z  axis  at  beam  i1-*1 
end. 

PLM35 

R*8 

150 

Plastic  moment  about  Y  axis  at  beam  i*-^1 
end. 

PLM26J 

R*8 

150 

Plastic  moment  about  Z  axis  at  beam 
end. 

PLM35J 

R*8 

150 

Plastic  moment  about  Y  axi.>  at  beam  j*1'1 
end . 

DEMA 

DEVMAX 

R*8 

Maximum  percent  deviation  for  mass  energy 
calculations. 
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TABLE  3-4.  COMMON  BLOCK  SUMMARY  (Continued) 


Common 

Name 


DINICP 


DOIN 

ENERGY 


Mnemonic 

Type 

Array 

Size 

XDOTAP 

R*8 

YDOTAP 

R*8 

ZDOTAP 

R*8 

DLVOC 

R*8 

5,  3 

FRD 

R*8 

130,4 

FRS 

R*8 

150,4 

SUMDFI 

R*8 

6, 150 

finti 

R*8 

6,150 

VEEN 

R*8 

2,150 

STENS 

R*8 

20 

SCOMP 

R*8 

20 

SHEAR 

R*8 

20 

EE 

R*8 

20 

GG 

R*8 

20 

FINT 

R*8 

6,150 

VOL 

R*8 

5 

VZERO 

R*8 

5 

KMATR 

1*4 

6,  4 

NVCH 

1*4 

INBUFF 

1*2 

Ln 

00 

CBAR 

R*8 

150 

XPCT 

R*8 

80 

XETOT 

R*8 

80 

XETOTO 

R*8 

80 

XSE 

R*8 

80 

XDE 

R*8 

80 

XCE 

R*8 

80 

XFE 

R*8 

80 

KEI 

R*8 

80 

Overall  vehicle  velocity 
parallel  to  ground  axis. 


Volume  length  changes. 

Ratio  of  current  max.  stress  to  failure 
stress  -  constant  energy  at  dist.  theory 

Ratio  of  current  max.  stress  to  failure 
stress  -  max.  shear  theory. 

Strain  force  at  i*-h  mass. 

Total  of  strain  and  damping  force  at  i1-*1 
mass . 

Sum  of  0  and  ip  deflections  at  ich  and  jtl 
ends  of  the  beam. 

Tension  yield  stress. 

Compression  yield  stress. 

Shear  yield  stress. 

Young's  modulus. 

Torsional  modulus. 

Total  beam  force. 


Current  "|  volume  for  volume  change 

Initial  J  calculations 

Volume  corner  numbers. 

Number  of  volumes. 

Mass  numbers  forming  volumes. 

Damping  coefficient  for  beam  IJ. 

Percent  of  total  energy  change  at  mass. 

Total  energy  (by  masses  and  for  system). 

Initial  total  energy  (by  masses  and  for 
system) . 

Mass  related  strain  energy. 

Mass  related  damping  energy. 

Mass  related  crushing  energy. 

Mass  related  friction  energy. 

Mass  related  kinetic  energy. 


TABLE  3-4.  COMMON  BLOCK  SUMMARY  (Continued) 


Common 

Name 

— 

Mnemonic 

Type 

Array 

Size 

ENERGY 

PEI 

R*8 

80 

Mass  related  potential  energy. 

(Cont' d) 

XETOTL 

R*8 

Total  system  energy. 

KETOTL 

R*8 

Total  system  kinetic 

energy . 

PETOTL 

R*8 

Total  system  potential  energy. 

SETOTL 

R*8 

Total  system  strain  energy. 

DETOTL 

R*8 

Total  system  damping 

energy . 

CETOTL 

R*8 

Total  system  crushing  energy. 

FETOTL 

R*8 

Total  system  friction  energy- 

IBALL 

IBUF1 

1*4 

20 

Index  for  new  materials. 

INAC 

ACCEL 

R*8 

300 

Table  of  option  mass 

acceleration  levels. 

TIM 

R*8 

300 

Time  at  which  acceleration  levels  are 

input . 

INDEX 

1*4 

50,  2 

Mass  number  and  direction  index  for 

acceleration  data. 

JAY 

1*4 

50,  2 

Table  number  and  number  of  acceleration 

points  in  table. 

KOUNT 

1*4 

Counter  for  number  of  acceleration  pulse 

points . 

INCF 

SA 

R*8 

40 

SB 

R*8 

40 

SF 

R*8 

40 

^  External  spring  input  table  parameters. 

SI 

R*8 

40 

d 

XMU 

R*8 

40 

External  spring  coeff.  of  friction. 

XKE 

R*8 

40 

Bottoming  spring  stiffeners. 

XMAX 

R*8 

40 

1 

plowing  force 

FSPOF 

R*8 

40 

-  External  spring  -1 

load-defl.  par. 

FSPOI 

R*8 

40 

J 

v  load-defl .  par . 

GFLEX 

R*8 

40 

Ground  stiffness 

CDAMP 

R*8 

40 

External  spring  damping  coefficient 

PLOW! 

R*8 

Starting  time  for  plowing  force. 

INCFIC 

BETA 

R*8 

Ground  slope  angle  - 

degrees. 

INDEAC 

NACC 

1*4 

Number  of  acceleration  input  tables. 

INIC 

XDP 

R*8 

80 

F 

■  S.  'I  ,  . 

L  coordinates 

jf  mass  points 

ZDP 

R*8 

80 

W.L.  J 
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TABLE  3-4.  COMMON  BLOCK  SUMMARY  (Continued) 


Common 

Array 

1 

Name 

Mnemonic 

Type 

Size 

INIC 

PHIDP 

R*8 

80 

N 

(Cont ' d) 

PSIDP 

R*8 

80 

Euler  angles  for 
axis  at  t  =  0.0. 

body  axis  to  airplane 

THEDP 

R*8 

80 

-> 

PPR 

R*8 

> 

QPR 

R*8 

”  Initial  airplane 

angular  velocities. 

RPR 

R*8 

- 

XGIN 

R*8 

Distance  aft  "1  of 

initial  c.g.  position 

ZGIN 

R*8 

Distance  up  j  of  vehicle 

PHIPR 

R*8 

' 

PSIPR 

R*8 

Initial  airplane 
to  ground. 

attitude  relative 

THEPR 

R*8 

y 

XGDOT 

R*8 

YGDOT 

R*8 

Initial  airplane 
^  ground  axis. 

c.g.  velocity  in 

ZGDOT 

R*8 

-■ 

INIDCP 

YDP 

R*8 

80 

Mass  "Y"  coordinate 

ININPR 

NSF 

1*4 

> 

*  beam  strain  forces 

NTF 

1*4 

total  beam  forces 

NDE 

1*4 

beam  deflections 

NSPD 

NED 

1*4 

1*4 

y.  Print  flags  for-^ 

spring  loads  and 
deflections 

NS 

1*4 

energy  breakdown  by 

mass 

NRP 

1*4 

beam  &  spring 

NIMP 

1*4 

J 

beam  stresses 

mass  displacements, 
velocities,  and 
accelerations 

s.  mass  impulses 

INOUT 

FCUT 

R*8 

Cutoff  frequency  for  digital  filter. 

NTOL1 

1*2 

Tolerances  for  total  energy  growth. 

NTOL2 

1*2 

Tolerance  for  negative  train,  damping, 
crushing  and  friction  length. 

NTOL3 

1*2 

Tolerance  for  mass  energy  deviation. 

INPR 

NDRI 

184 

Number  of  DRI  beam  elements. 

NSP 

1*4 

Number  of  external  springs. 
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TABLE  3-4.  COMMON  BLOCK  SUMMARY  (Continued) 


Common 

Array 

Name 

Mnemonic 

Type 

Size 

IPIC 

NIC 

1*4 

Flag  to  call  for  print  out  of  preliminary 
failure  loads  and  deflections. 

IRDE 

PFIL 

R*8 

Tire  constant. 

MACF 

FSOP 

R*8 

40 

Maximum  spring  force  -  previous  value. 

SCP 

R*8 

40 

Maximum  spring  deflection  -  previous 
value . 

1.  Initial  deflection  value. 

2.  Maximum  force  value. 

EXSP 

R*8 

5,200 

4 

3.  Maximum  deflection  value. 

4.  Unloaded  deflection  value. 

► 

5.  Unloaded  force  value. 

^  J 

TSP 

R*8 

200 

r  Time  for  external  spring  loading 
.  or  unloading  event. 

STEMP1 

R*8 

40 

STEMP2 

R*8 

40 

*< 

Temporary  tracking  of  external  spring 

- 

STEMP3 

R*8 

40 

deflections  and  forces  for  use  in 

STEMP4 

R*8 

40 

restart . 

STEMP5 

R*8 

40 

r 

1.  Spring  mass  number. 

2.  Massless  node  number. 

I 

KKSP 

1*4 

5,  200 

< 

3.  Spring  direction. 

4.  Type  =  JBS. 

► 

5.  Not  used. 

1 

KKONT 

1*4 

Total  number  of  liaes  of  print  for  exter¬ 
nal  spring  behavior  summary. 

IBS 

1*2 

40 

Flag  cut  to  1  when  external  spring  has 
bottomed. 

JBS 

1*2 

40 

Flag  used  in  summary  of  external  spring 
behavior : 

1.  Initial  loading. 

2.  Max.  loading. 

3.  Unload  to  zero. 

4.  Start  of  reload. 

MACFIN 

IPRINT 

1*4 

Multiple  of  time  increment  defining  print 

time . 

TMAX 

R*8 

Maximum  run  time. 

MADE 

KFL26 

1*4 

150 

Flag  that  indicates  whether  direction  2  or 

f 

went  nonlinear  first. 

KFL35 

1*4 

150 

Flag  that  indicates  whether  direction  3  or 

went  nonlinear  first. 
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TABLE  3-4.  COMMON  BLOCK  SUMMARY  (Continued) 


Common 

Name 

MADE 

(Cont' d) 


Mnemonic  Type 


FL26I  1*2 
FL26J  1*2 
FL35I  1*2 
FL35J  1*2 


ETOTTO  I R*8 


Array 

Size 


100,3 


Flag  to  indicate  occurrence  of  plastic 
hinge  formation  about  Z  axis  at  i*^1  & 
jth  end  and  about  Y  axis  at  i1-*1  and 
j th  end,  respectively. 


Plastic  hinge  summary  beam  number 
direction,  and  ic^  or  jth  end. 

Counter  for  number  of  plastic  hinge 
formations . 

Time  of  plastic  hinge  formation. 
Total  system  energy  for  time  zero. 


ENGSMY  R*4  6,200  Energy  summary  terms  printed  at  end  of 


1APR1 


INGSCT  1*4 


ETIME 


IEPSE 


IEPDE 


IEPCE 


IEPFE 

IEDEV 

IETOT 

IRUPSW 

IPENSW 


VEEBAR  R*8 
ZINIT  R*8 


DPHIIJ  R*8 


Total  number  of  time  cuts  printed  in 
energy  summary. 

Absolute  value  of  mass  energy  deviation. 
Summary  table  for  energy  print. 

Error  flag. 

x  f'total  strain  energy 

element  strain  energy 

total  damping  energy 

Error  element  damping  energy 

J  total  crushing  energy 
V  flags  -< 

^  mass  crushing  energy 

total  friction  energy 
mass  friction 
mass  deviation 

✓  l^total  energy  change 

Identification  of  "rupture"  beams. 

Identifies  mass  volume  penetration 
occurrence . 

Absolute  value  of  beam  deflection. 

Initial  value  (time  =  0)  of  mass  ground 
"Z"  coordinate  value. 

Incremental  value  of  Euler  angle  for  $ 
direction .' 
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TABLE  3-4.  COMMON  BLOCK  SUMMARY  (Continued) 


Common 

Array 

Name 

Mnemonic 

Type 

Size 

MAX 

FUB 

R*8 

150 

Total  axial  load 

for  an  unsymmetr ical 

(Cont ' d) 

beam. 

MCFIII 

SYMFLG 

R*8 

Flag  defining  symmetry  option. 

NPOOI2 

MG 

1*2 

150 

Node  point  number  at  i*-^1  end  of  beam. 

NG 

1*2 

150 

Node  point  number  at  j end  of  beam. 

INP 

1*2 

50 

Mass  to  which  node  is  rigidly  attached. 

MNP 

1*2 

50 

Node  point  number. 

NPOOI4 

NMP 

1*4 

Number  of  mass  less  node  elements. 

NPOOR8 

RX 

R*8 

50 

X1 

RY 

R*8 

50 

Y 

component  of  vector  from  mass  to  node 
point,  mass  axes 

RZ 

R*8 

50 

zj 

NP01I1 

II 

1*2 

40 

Mass  number  for  ^ 

KK 

1*2 

40 

Direction  of 

«■  IKMtl1  spring 

MM 

1*2 

40 

Node  point  for  , 

NP01R8 

XNP 

R*8 

50 

X' 

YNP 

R*8 

50 

Y 

►coordinates  of  node  points 

ZNP 

R*8 

50 

zj 

UNP 

R*8 

50 

X1 

VNP 

R*8 

50 

component  of  node  point  velocity, 
mass  axes 

MNP 

R*8 

50 

zj 

XNDP 

R*8 

50 

X1 

YNDP 

R*8 

50 

Y 

component  of  node  point  velocity, 
ground  axes 

ZNDP 

R*8 

50 

zj 

XACCNP 

R*8 

50 

XI 

YACCNP 

R*8 

50 

vl 

component  of  node  point  acceleration, 
mass  axes 

ZACCNP 

R*8 

50 

zj 

SBUCKR 

R*8 

150 

Beam  critical  buckling  stress 

PCR 

R*8 

150 

Beam  critical  buckling  load 

NP02I2 

MQ 

1*2 

180 

\ 

Node  points  for  spec,  of  KR  tables. 

NQ 

1*2 

180 

J 

IJVM 

1*2 

150 

1 

Index  for  beams  which  have  positive 

[ 

deflection  or 

force  for  a  rupture 

IJFM 

1*2 

150 

J 

value,  respectively. 
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TABLE  3-4.  COMMON  BLOCK  SUMMARY  (Continued) 


Common 

Name 

Mnemonic 

Type 

Array 

Size 

NP02I2 

IJVMN 

1*2 

150 

"3 Index  for  beams  which  have  negative 

Vdeflection  or  force  for  a  rupture 

value , 

IJFMN 

1*2 

150 

J  respectively. 

NP02R8 

XNPDP 

R*8 

50 

-\ 

r  F.S. 

YNPDP 

R*8 

50 

>-  Node  point  coordinates 

B.L. 

ZNPDP 

R*8 

50 

> 

„  W.L. 

OLEO 

EOLEO 

R*8 

20 

Effective  total  strut  cylinder 

le 

ngth. 

FAO 

R*8 

20 

Fully  extended 

gear  preload 

FAA 

R*8 

20 

Ambient  air  preload. 

EXPOLE 

R*8 

20 

Polytropic  exponent 

YMAX 

R*8 

20 

Maximum  strut  stroke 

YOLEO 

R*8 

20 

Strut  stroke  at  any  time. 

BOLEO 

R*8 

20 

Strut  orifice  damping. 

BROLEO 

R*8 

20 

Strut  rebound  valve  damping 

XKEXT 

R*8 

20 

Linear  spring  stiffness  at  extended  end 

of  strut  travel. 

XKCOMP 

R*8 

20 

Linear  spring  stiffness  at  compressed  end 

of  strut  travel. 

FCOUL 

R*8 

20 

Coulomb  or  constant  friction  force  value. 

ALPHA? 

R*8 

Constant  used  in  computing  shock 

strut 

friction  force 

IGOLEO 

1*4 

20 

i 

th  end  of  beam  identification 

• 

JGOLEO 

1*4 

20 

j 

th  end  of  beam  identification 

MGOLEO 

1*4 

20 

m 

massless  node  identification. 

NGOLEO 

1*4 

20 

n 

th  massless  node  identification. 

NOLEO 

1*4 

Number  of  shock  struts. 

OTPLT 

NMEP 

1*2 

'’mass  > 

NNEP 

1*2 

> 

node 

NBFP 

1*2 

beam  force 

NBDP 

1*2 

beam  deflec- 

plots 

K  Number  of  A 

y 

NSEP 

tion  spring 

requeste 

NDRP 

1*2 

DRI 

NSTP 

1*2 

beam  stress 

NENP 

1*2 

J 

v beam  energy  ^ 
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TABLE  3-4.  COMMON  BLOCK  SUMMARY  (Continued) 


Common 

Name 

Mnemonic 

Type 

Array 

Size 

i -  - 1 

i _  1 

OTPLT 

JMASS 

1*2 

50, 10 

f  mass 

\ 

(Cont’ d) 

JNODE 

1*2 

50,8 

node 

JBMF 

1*2 

50,4 

beam  force 

JBMD 

1*2 

50,4 

Plot 

beam  deflec- 

>-  flag 

t  ion 

L  plot 

JBMS 

1*2 

50,6 

for 

beam  stress 

JSPR 

1*2 

50,4 

spring 

JENG 

1*2 

50,3 

beam  energy 

JDRI 

1*2 

10 

J 

c  DRI 

J 

NMEW 

1*2 

> 

'"mass 

NNEW 

1*2 

node 

NBFW 

1*2 

beam  force 

Counters  which  are 

.  . . 

NBDW 

1*2 

L 

beam  deflection 

\  indexed  when  plot  ^ 

NStW 

1*2 

spring 

data  is  written 

NDRW 

1*2 

DRI 

NSTW 

1*2 

beam  stress 

NENW 

1*2 

^beam  energy 

NPRINT 

1*2 

Flag 

which  controls  counters. 

PPLTS 

XSCALE 

R*8 

10 

X 

axis  f 

scale  factor  for 

Y SCALE 

R*8 

10 

Y 

axis  J 

position  plots 

NPLT 

1*2 

Number  of  plots. 

NPFCT 

1*2 

Plot 

print  factor. 

IPFCT 

1*2 

Counter  for  plot  print. 

ITPL 

1*2 

10 

Plot 

plane . 

NMPTS 

1*2 

10 

Number  of  points  in  plot. 

MNUM 

1*2 

50, 10 

Mass 

point  identification  of  plot  points. 

ISCALE 

1*2 

10 

Scale  option. 

|  PRMA 

XACF 

R*8 

80 

r 

YACF 

R*8 

80 

-  ■< 

Mass  point  filtered  acceleration. 

ZACF 

R*8 

80 

> 

L 

XACNPF 

R*8 

50 

I 

r 

YACNPF 

R*8 

50 

i 

Node  point  filtered  acceleration. 

I _ 

ZACNPF 

R*8 

50 

J 

1 
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TABLE  3-4 .  COMMON  BLOCK  SUMMARY  (Continued) 


Common 

Array 

— 

Name 

Mnemonic 

Type 

Size 

PRMA 

XIMP 

R*8 

80 

r 

(Cont ' d) 

YIMP 

R*8 

80 

1 

Mass  point  filtered  impulses. 

ZIMP 

R*8 

80 

L 

XIMPNP 

R*8 

50 

r 

YIMPNP 

R*8 

50 

- 

Massless  node  point  filtered  impulses. 

ZIMPNP 

R*8 

50 

XIMPOL 

R*8 

80 

r 

YIMPOL 

R*8 

80 

- 

Mass  point  filtered  impulses  - 
previous  value 

ZIMPOL 

R*8 

80 

i 

XIMPNL 

R*8 

50 

r 

! 

YIMNPL 

R*8 

50 

Massless  node  point  filtered  impulses  - 
previous  value. 

ZIMNPL 

R*8 

50 

i 

RESTRT 

CASEIN 

R*8 

Checkpoint  title. 

RUN  IN 

1*4 

Check  point  case  number. 

MSEC IN 

1*4 

Check  point  save  times. 

CASOUT 

R*8 

Restart  title. 

RUNOUT 

1*4 

Restart  case  number. 

MSCOUT 

1*4 

5 

Restart  time. 

STUFF 

NEWI 

1*4 

80 

f 

NEWIJ 

1*4 

150 

►  < 

Beam  and  mass  numbering  correspon- 

NMOLD 

1*4 

dence  (left/right)  data. 

NBOLD 

1*4 

- 

UB 

DB 

R*8 

150 

Deadband. 

IJUB 

1*4 

150 

Type; 

1  -  tension 

-1  -  compression 

NUB 

1*4 

Number  of  unsymmetrical  beams. 

**VARINT 

MINDT 

R*8 

Minimum  integration  interval. 

DT2 

R*8 

Half 

of  integral  interval. 

TPRINT 

R*8 

Print  time. 

EL 

R*8 

Lower  limit. 

EU 

R*8 

Upper  limit. 

i _ 

RATMIN 

R*8 

Minimum  rate  of  change  of  DT. 
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3.3  PERMANENT  DATA  STORAGE  AND  RETRIEVAL 


The  tape  data  management  routines  included  in  program  KRASH  are  executed 
when  the  restart  option  is  invoked.  The  routines  provide  a  systematic  means 
by  which  data  can  be  stored  on  and  retrieved  from  magnetic  tape.  A  general 
description  of  the  management  system  and  routines  involved  is  provided  in  the 
following  sections. 

3.3.1  Tape  Data  Management 

3.3. 1.1  Data  Storage 

Before  any  data  can  be  stored,  it  is  necessary  to  first  place  an  end  of 
file  (EOF)  mark  at  the  beginning  of  the  reserved  tape  allocated  for  the 
permanent  storage  of  restart  data.  This  is  necessary  because  when  data  are 
stored  the  routines  search  for  an  EOF  mark  before  storing. 

Data  are.stored  in  blocks  which  are  referred  to  as  "user  data  sets". 

These  data  sets  are  made  up  of  the  state  vector  data  referred  to  in  Section  2.3. 
Prior  to  storing  a  user  data  set,  a  header  is  generated  and  placed  at  the 
beginning  of  the  user  data  set  along  with  restart  information  supplied  as 
input  which  uniquely  identifies  the  data  set.  See  Section  2.1  of  Reference  (2) 
and  line  6  of  the  Example  Input  Data,  Figure  4-2  of  this  manual. 

3. 3. 1.2  Data  Retrieval 

A  user  data  set  previously  stored  as  part  of  the  restart  option  may  be 
retrieved  by  searching  the  stored  data  for  the  header  information  which 
identifies  the  desired  data  set.  The  header  information  identifying  the 
desired  data  set  is  supplied  as  input  data.  See  Section  2.1  of  Reference  (2) 
and  line  5  of  the  Example  Input  Data,  Figure  4-3  of  this  manual. 

3.3.2  Data  Management  Routines 

3. 3. 2.1  Data  Storage 

3. 3.2. 1.1  OPNOUT 


OPNOUT  is  an  entry  point  in  subroutine  DATOUT. 
made  from  subroutine  RSOUT  to  open  a  new  data  set. 


The  call  to  OPNOUT  is 


OPNOUT  reads  through  the  data  on  the  restart  tape  until  an  EOF  mark  is 
located.  The  tape  is  backspaced  over  the  EOF  mark  and  the  user  data  set  mark 
(**START*)  and  header  is  written  over  the  EOF  mark.  The  header  consists  of  a 
24-word  (96  byte)  record  containing  the  following  information: 

•  Words  1-2  Date  in  the  form  MM-DD-YY  supplied  by  a  call  to 
subroutine  CDATE 


•  Words  3-24  Data  supplied  by  a  call  to  subroutine  JOBAC 


CALLING  STATEMENT  ARGUMENT  LIST 

Mnemonic 

— ,  - 

Symbol 
R:-f.  (1) 

Type 

UNIT 

n 

Unit  allocated  for  the  reversed  restart 
tape  (Unit  21) 

3.3.2. 1.2  DATOUT 

Subroutine  DATOUT  is  called  by  subroutine  RSOUT .  Data  are  collected 
within  DATEOUT.  To  these  data  are  added  an  8-character  literal  used  to  name 
the  data  set  and  a  four  byte  word  which  specifies  the  number  of  words  to  be 
stored.  Thus,  the  total  data  set  consists  of  (1)  an  8-character  literal, 

(2)  a  4  byte  word,  and  (3)  the  state  vector  data.  All  data  are  stored  in 
blocks  of  250  words  in  core.  When  the  block  is  filled,  the  data  are  written 
onto  the  reserved  restart  tape. 


CALLING  STATEMENT  ARGUMENT  LIST 


Mnemonic 

Symbol 
Ref.  (1) 

Type 

X 

State  vector  data  to  be  saved 

NAME 

NAME 

Check  point  run  name  identifier 

NWDS 

Number  of  words  to  be  stored 

UNIT 

Unit  allocated  for  the  reserved  restart 

tape  (Unit  21) 

ERROR 

Error  flag  returned  by  the  tape  write 

mm 

routine 
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3.3.2. 1.3  CLSOUT 


CLSOUT  is  an  entry  point  in  subroutine  DATOUT.  The  call  to  CLSOUT  is 
made  from  subroutine  RSOUT  to  close  a  new  data  set. 

CLSOUT  fills  any  remaining  space  in  the  250-word  block  used  by  DATOUT 
with  end  of  data  set  marks  (*END) .  After  the  block  is  filled  and  written  onto 
the  reserved  restart  tape  an  EOF  mark  is  placed  on  the  tape.  The  tape  is 
then  backspaced  over  the  EOF  mark  in  anticipation  of  a  further  extension  of 
the  tape. 


CALLING  STATEMENT  ARGUMENT  LIST 

Mnemonic 

Symbol 
Ref.  (1) 

Type 

UNIT 

1*4 

Unit  allocated  for  the  reserved  restart 
tape  (Unit  21) 

3. 3. 2. 1.4  JOBAC 

Subroutine  JOBAC  is  supplied  as  a  dummy  subroutine.  It  provides  a  means 
by  which  user  supplied  accounting  or  other  pertinent  information  can  be 
returned  to  the  calling  subroutine  DATOUT  as  a  22-word  array  named  ACCARD. 

As  supplied,  JOBAC  returns  blank  data. 


CALLING  STATEMENT  ARGUMENT  LIST 

Mnemonic 

Symbol 
Ref.  (1) 

Type 

ACCARD 

R*4 

Twenty-two  word  blank  array 
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3.3.2. 1.5  CDATE 


Subroutine  CDATE  is  supplied  as  a  dummy  subroutine.  It  provides  a  means 
by  which  the  current  date,  to  be  supplied  by  the  user,  can  be  returned  to 
the  calling  subroutine  DATOUT  by  way  of  the  calling  statement  agreement  list. 
As  supplied,  this  subroutine  returns  the  character  strings  MM-D  and  D-YY. 


CALLING  STATEMENT  ARGUMENT  LIST 

Mnemonic 

Symbol 
Ref.  (1) 

Type 

DATE 

R*4 

Dummy  array  name 

ITYPE 

1*4 

Dummy  name 

3. 3. 2. 2  Data  Retrieval 
3. 3. 2. 2.1  OPNIN 

OPNIN  is  an  entry  point  in  subroutine  DATIN.  The  call  to  OPNIN  is  made 
from  subroutine  RSIN  to  open  an  existing  user  data  set  prior  to  the  retrieval 

of  data. 


OPNIN  reads  the  reserved  restart  tape  until  a  user  data  set  mark  (**START*) 
is  found.  The  first  250  word  block  after  the  user  data  set  mark  is  read  and 
the  header  and  data  set  identifiers  are  recovered  and  returned  to  RSIN. 


CALLING  STATEMENT  ARGUMENT  LIST 

Symbol 

Mnemonic 

Ref.  (1) 

Type 

HEADER 

1*4 

Header  array 

NWDS 

Number  of  words  to  be  retrieved 

UNIT 

Unit  allocated  for  the  reserved  restart 

■Ml 

tape  (Unit  21) 

ERROR 

n 

Error  flag  returned  by  tape  read  routine 
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3. 3. 2. 2. 2  DATIN 


Subroutine  DATIN  Is  called  from  subroutine  RSIN.  DATIN  receives  header 
and  user  data  set  identifiers  from  OPNIN.  The  identifiers  are  compared  with 
corresponding  identifiers  supplied  as  input  data-  See  Section  2.1  of 
Reference  (2).  If  agreement  is  achieved,  the  remaining  data  in  the  user  data 
set  are  read  from  the  reserved  restart  tape  and  returned  to  RSIN.  If  the 
comparison  fails,  execution  is  passed  to  CLSIN  where  the  data  set  is  closed 
and  OPNIN  called  again. 


CALLING  STATEMENT  ARGUMENT  LIST 

Mnemonic 

Symbol 
Ref.  (I) 

Type 

X 

R*8 

State  vector  quantity  retrieved 

NAME 

R*8 

Restart  run  name  identifier 

NWDS 

1*4 

Number  of  words  to  be  returned 

UNIT 

1*4 

Unit  allocated  for  the  reserved  restart 
tape  (Unit  21) 

ERROR 

1*4 

Error  flag  returned  by  tape  read 
rout ine 

3. 3. 2. 2. 3  CLSIN 


CLSIN  is  an  entry  point  in  subroutine  DATIN.  The  call  to  CLSIN  is  made 
from  subroutine  RSIN  to  close  a  currently  opened  user  data  set  by  scanning  the 
data  set  until  an  end  of  data  mark  (*END)  is  found. 


CALLING  STATEMENT  ARGUMENT  LIST 

Mnemonic 

Symbol 
Ref.  (1) 

Type 

UNIT 

1*4 

Unit  allocated  for  reserved  restart 

ERROR 

B 

tape 

Error  flag  returned  by  tape  read  routine 
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3. A  ERROR/WARNING  MESSAGES 


Tests  of  pertinent  data  are  made  at  strategic  locations  throughout  progra 
KRASH.  These  checks  are  designed  to  identify  erroneous  input  data  and  incon¬ 
sistent  system  operations  and  computational  results.  A  summary  of  the  error 
messages  resulting  from  the  failure  of  any  of  these  tests  is  presented  in 
Table  3-5.  Referenced  in  the  table  is  the  subroutine  and  format  statement 
generating  the  message. 


1 


TABLE  3-5.  ERROR/WARNING  MESSAGE  SUMMARY 


Message 

No. 

Error/Warning  Message 

Format 

Label 

Sub- 

Routine 

1 

IDER  OUT  OF  RANGE,  VALUE  = 

Degree  of  freedom  flag  assigned  to  tabu¬ 
lated  input  acceleration  data  not  within 
the  allowable  range  of  1  through  6. 

Control  returned  to  calling  routine. 

6106 

ACCELT 

2* 

NO  SUCH  NODE  POINT  I,  M  = 

External  spring  attached  to  node  which 
has  not  previously  been  defined.  Run 
terminated . 

9 

CFORCE 

3* 

DATA  ERROR  IN  GENMOD  AT  2030 

Unable  to  find  specified  node  point.  Run 
terminated . 

2030 

GENMOD 

4* 

DATA  ERROR  IN  GENMOD  AT  2040 

Unable  to  find  specified  node  point.  Run 
terminated . 

2040 

GENMOD 

5* 

NON-EXISTENT  I,  J  PAIR  IN  KR  TABLE 

SPECS  = 

Unable  to  locate  beam  I,  J  in  KR  table. 

Run  terminated. 

5910 

INPUT 

6* 

TOO  MANY  KR  TABLES 

Number  of  KR  tables  exceeds  program 
limit.  Run  terminated. 

5920 

INPUT 

7 

ERROR  IN  ACCEL  TIME  DATA  INPUT  NUMACL, 

KOUNT  = 

Number  of  tabulated  acceleration  tables 
or  the  total  number  of  tabulated  accel¬ 
eration  points  exceeds  the  program 
limit.  Run  execution  continues. 

7005 

INPUT 

8* 

POINTS  IN  KR  TABLE  (MAX  IS  15) 

Too  many  points  have  been  specified  in 

KR  table.  Run  terminated. 

5980 

INPRNT 

9* 

NUMBER  OF  LINEAR  DEFLECTION  POINTS  NOT  AC¬ 
CEPTABLE,  VALUE  - 

Number  of  linear  deflection  points  (LDP) 
defining  KR  curve  exceeds  program  limit. 

Run  terminated. 

5088 

INPRNT 

10 

I/O  ERROR  ON  READ  -  WILL  REWIND  INPUT  FILE 

AND  RETURN. 

Run  continues. 

60 

ECHO 

Error  causing  run  to  terminate  execution 
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TABLE  3-5.  ERROR/WARNING  MESSAGE  SUMMARY  (Continued) 


Message 

No. 

Error/Warning  Message 

Format 

Label 

Sub- 

Routine 

11* 

****ERROR  -  END  OF  DATA  REACHED  AT 

VARIABLE  **** 

A  greater  amount  of  data  requested  from 
restart  tape  than  available.  Run 
terminated . 

1000 

DAT  IN 

12* 

***ERROR  -  END  OF  FILE  REACHED  AT  VARIABLE 

— — 

The  data  set  currently  read  from  restart 
tape  was  not  closed  after  it  was  created. 
Run  terminated. 

3000 

DATIN’ 

13* 

****ERROR  -  END  OF  FILE  REACHED  TRYING  TO 
OPEN  DATA  SET  FOR  READING 

The  restart  tape  does  not  contain  a 
previously  saved  data  set.  Run 
terminated . 

4000 

DATIN' 

14* 

****ERROR  -  EOF  TRYING  TO  CLOSE  DATA  SET 
AFTER  READING**** 

The  data  set  currently  read  from  restart 
tape  was  not  closed  after  it  was  created. 
Run  terminated. 

5000 

DATIN 

15* 

****ERROR  -  CANNOT  READ  DATA  WITHOUT  FIRST 
OPENING  DESIRED  DATA  SET**** 

The  data  set  to  be  read  from  restart 
tape  was  not  opened  prior  to  request. 

Run  terminated. 

6000 

DATIN 

16* 

****ERROR  -  MUST  OPEN  DATA  SET  BEFORE  WRIT¬ 
ING  DATA  OUT**** 

Data  set  not  opened  prior  to  writing  on 
restart  tape.  Run  terminated. 

1000 

DATOUT 

17 

**WARNTNG  -  REOUESTED  STORING  VALUES. 
STORE  FAILED** 

The  number  of  words  to  be  written  on  the 
restart  tape  exceeds  the  space  allocated 
by  the  program.  Run  continues. 

2000 

DATOUT 

18 

UNABLE  TO  FIND  DATA  FOR  MASS  .  PLOTS 

FOR  THIS  MASS  WILL  BE  SUPPRESSED. 

Data  for  indicated  mass  summary  plots  not 
found  on  unit  1.  Run  continues. 

6023 

PREPLT 
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TABLE  3-5.  ERROR/WARNING  MESSAGE  SUMMARY  (Continued) 


Message 

No. 

Error /Warning  Message 

Format 

Label 

Sub- 

Routine 

19 

UNABLE  TO  FIND  DATA  FOR  MASS  u(W 

.  PLOTS  FOR  THIS  NODF.  WILL  RE 

SUPPRESSED. 

Data  for  indicated  node  point  summary 
plots  not  found  on  unit  2.  Run 
continues. 

6024 

PREPLT 

20 

UNABLE  TO  FIND  DATA  FOR  BEAM  .  FORCE 

PLOTS  FOR  THIS  BEAM  WILL  BE  SUPPRESSED. 

Data  for  indicated  beam  summary  force 
plots  not  found  on  unit  3.  Run  continues. 

6025 

PREPLT 

21 

UNABLE  TO  FIND  DATA  FOR  BEAM  .  DEFT. EG- 

TION  PLOTS  FOR  THIS  BEAM  WILL  BE  SUPPRESSED. 
Data  for  indicated  beam  summary  deflection 
plots  not  found  on  unit  4.  Run  continues. 

6026 

PREPLT 

22 

UNABLE  TO  FIND  DATA  FOR  BEAM  .  STRESS 

PLOTS  FOR  THIS  BEAM  WILL  BE  SUPPRESSED. 

Data  for  indicated  beam  summary  stress 
plots  not  found  on  unit  13.  Run 
continues . 

6029 

PREPLT 

23 

UNABLE  TO  FIND  DATA  FOR  SPRTNG  NODE 

_  PLOTS  FOR  THIS  SPRING  WILL  BE 

SUPPRESSED. 

Data  for  indicated  spring  summary  plots 
not  found  on  unit  8.  Run  continues. 

6027 

PREPLT 

24 

UNABLE  TO  FIND  DATA  FOR  DRI  MASS 

PLOTS  FOR  THIS  DRI  MASS  WILL  BE  SUPPRESSED. 
Date  for  indicated  DRI  mass  summary  plots 
not  found  on  unit  9.  Run  continues. 

6028 

PREPLT 

25 

***F.RROR  PLOT***  T  «=  _  NOT  « 

Error  in  reading  vehicle  velocity  data 
from  unit  10.  Run  continues. 

1000 

PREPLT 

26 

UNABLE  TO  FIND  DATA  FOR  BEAM  .  ENERGY 

PLOTS  FOR  THIS  BEAM  WILL  BE  SUPPRESSED. 

Data  for  indicated  beam  energy  plots  not 
found  on  unit  11.  Run  continues 

6040 

PREPLT 

27* 

IESE  =  1  TOTAL  STRAIN  ENERGY  IS  NEGATIVE. 

The  total  strain  energy  associated  with 
the  model  is  less  than  aero. 

7301 

MAIN 
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TABLE  3-5.  ERROR/WARNING  MESSAGE  SUMMARY  (Continued) 


Message 

No. 

Error /Warning  Message 

Format 

Label 

28* 

IEPSE  =  1  ELEMENT  STRAIN  ENERGY  EXCEEDS 

NT0L2  VALUE. 

The  strain  energy  of  any  element  exceeds 
a  prescribed  percentage  of  the  total 
strain  energy  of  the  model. 

7302 

29* 

IEDE  =  TOTAL  DAMPING  ENERGY  IS  NEGATIVE. 

The  total  damping  energy  associated  with 
the  model  is  less  than  zero. 

7303 

30* 

IEPDE  =  1  ELEMENT  DAMPING  ENERGY  EXCEEDS 

NT0L2  VALUE. 

The  damping  energy  of  any  element  ex¬ 
ceeds  a  prescribed  percentage  of  the 
total  damping  energy  of  the  model. 

7304 

31* 

IECE  =  1  TOTAL  CRUSHING  ENERGY  IS  NEGATIVE. 

The  total  crushing  energy  associated  with 
the  model  is  less  than  zero. 

7305 

32* 

IEPCE  =  1  MASS  CRUSHING  ENERGY  EXCEEDS  NT0L2 

VALUE. 

The  crushing  energy  associated  with  any 
mass  exceeds  a  prescribed  percentage  of 
the  total  crushing  energy  of  the  model. 

7306 

33* 

IEFE  =  1  TOTAL  FRICTION  ENERGY  IS  NEGATIVE. 

The  total  friction  energy  associated  with 
the  model  is  less  than  zero. 

7307 

34* 

IEPFE  =  1  MASS  FRICTION  ENERGY  EXCEEDS  NT0L2 

VALUE. 

The  friction  energy  associated  with  any 
mass  exceeds  a  prescribed  percentage  of 
the  total  friction  energy  of  the  model. 

7308 

35* 

IEDEV  =  1  MASS  DEVIATION  EXCEEDS  NT0L3 

VALUE . 

The  total  energy  associated  with  any  mass 
exceeds  the  total  initial  energy  asso¬ 
ciated  with  that  mass  by  a  prescribed 
percentage . 

7309 

36* 

IETOT  -  1  TOTAL  ENERGY  CHANGE  EXCEEDS  NT0L1 
VALUE. 

The  total  energy  associated  with  the 
model  exceeds  the  total  initial  energy  by 
a  prescribed  percentage. 

7310 

Sub- 

Routine 

MAIN 


MAIN 

MAIN 


MAIN 

MAIN 


MAIN 

MAIN 


MAIN 


MAIN 
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SECTION  4 


DEMONSTRATION  PROBLEM 


The  example  problem  included  in  this  section  has  been  designed  to 
meet  two  objectives.  The  first  is  to  illustrate  an  input  data  card  deck 
which  exercises  some  commonly  used  options  available  to  the  user  of  the 
program.  The  second  is  to  provide  output  data  to  be  used  as  an  operational 
and  numerical  check  of  the  program.  The  model  used  in  the  example  does 
not  represent  any  particular  vehicle  nor  any  particular  crash  condition. 

4.1  MODEL  DESCRIPTION 

The  left  side  of  the  model  developed  for  this  example  is  shown  in 
Figure  4-1.  The  composite  model  consists  of  20  external  springs  (NSP), 

4  nonlinear  beams  (NLB),  12  massless  node  points  (NNP),  10  pin-ended 
beams  (NPIN),  8  unsymmetrical  beams  (NUB),  2  DRI's  (DRl) ,  2  beam  damping 
inputs  (ND)  and  a  RUNMOD  =2,  which  is  a  half  model  input  and  subsequent 
full  vehicle  analysis.  The  impact  is  with  a  flat  surface  at  a  vertical 
sink  velocity  of  27.5  ft/sec  and  a  0.75“degree  nose-up  pitch  angle.  The 
computer  run  is  performed  with  a  1  x  10  ^  integration  interval  and  a 
maximum  time  of  2  x  10  ^  seconds.  The  RUNMOD  =  2  is  chosen  to  demonstrate 
how  centerline  and  opposite  side  members  are  automatically  selected 
by  program  KRASH.  The  impact  condition  is  symmetrical,  so  for  this 
type  of  analysis  the  user  will  normally  use  the  RUNMOD  =  1  option,  which 
is  more  economical  to  perform. 

4.2  INPUT  DATA  DECK 

The  input  data  deck  used  in  the  demonstration  run  is  shown  in 
Figure  4-2.  The  card  identification  numbers  (columns  73-80)  used  in 
this  deck  does  not  follow  the  numbering  convention  used  in  the  program 
KRASH  User's  Manual,  Figure  2-1  of  Reference  (2).  Use  of  the  system 
shown  in  the  User's  Manual  is  not  mandatory.  Note  that  the  last  card 
of  the  input  deck  contains  the  work  "END"  beginning  in  column  one. 
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Figure  4-3  Example  Input  Data-Checkpoint  and  Restart  Runs 


For  illustrative  purposes,  Figure  4-3  shows  a  portion  of  the  UCHt 
decks  for  a  checkpoint  run  which  generates  data  for  a  future  resta, • 
run  and  a  restart  run  which  initiates  at  the  previsouly  designated 
save  time.  Output  from  these  latter  two  data  cases  are  not  included 
in  Appendix  A. 

4.3  OUTPUT 

A  complete  listing  of  the  output  for  the  demonstration  run  is 
included  in  this  report  as  Appendix  A.  The  format  of  the  data  is 
discussed  in  detail  in  Section  2.2  of  the  User's  Manual. 
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